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Aspice convexo nutantem pondere mundum 
terrasque tractusque maris caelumque profundum. 

Behold the world swaying her convex mass, 
lands and spaces of sea and depth of sky. 

(Vergilius, Ecloga IV^) 



^Translated from Latin by J W MacKail, in Virgils' Works (Modern Library, New York, 
1934). 
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Preface Note 



This manuscript must be intended as an informal review of the research works 
carried out during three years of PhD. "Informal" in the sense that technical 
proofs are often omitted (they can be found in the papers) as one could do for a 
presentation in a public talk. Clearly, some background of Quantum Mechanics is 
needed, even if I tried to minimize the prerequisites. 
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Introduction 



To handle information is to handle physical systems, and viceversa. Hence, the 
ultimate limits in manipulating and distributing information are posed by the very 
laws of physics. This is true in the classical framework (e. g. Landauer's principle) 
and in the quantum framework, where the rules of Quantum Mechanics give rise 
to new — and often not yet completely understood — restrictions and advantages 
to information processing and distribution. Quantum Information Theory is de- 
voted to the investigation of the theoretical limits Quantum Mechanics establishes 
when dealing with information encoded on quantum systems. This thesis treats 
the problem of processing quantum information^ in an optimal way by means of 
physically realizable devices. In fact, linearity of Quantum Mechanics forbids ba- 
sic processings of classical information — like e. g. copying-, broadcasting-, and 
NOT-gates — to properly work on an unknown quantum state. The first natural 
question is then: How well can we approximate such transformations and which 
are the physical devices that realize these approximations? 

In contrast to its classical counterpart, quantum information is very sensitive 
to noise. In a realistic setup, it is unreasonable to completely rule out noise, since 
the least interaction with the sorroundings can cause the system to be irreversibly 
disturbed. This fact raises the need of designing methods to encode quantum 
information in a way that is robust with respect to noise. But in order to do this, 
we have to provide a model for the noise. In this sense, also noise can be viewed as 
a kind of processing of quantum states: a "nasty" processing, nonetheless obeying 
the same laws of Quantum Mechanics as "good" processings do. The second 
natural question is then: What is the role of noise in a realistic setup and how can 

^Here "quantum information" is a short hand for "information encoded on quantum systems" 
and it is basically equivalent to saying "quantum states". Analogously, "classical information" 
means "information encoded on classical systems" . 
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we control it? 

In order to answer both questions, we clearly need to work in full general- 
ity. The appropriate mathematical tool to do this is provided by the concept of 
quantum channel. It encloses all possible deterministic transformations of quan- 
tum states allowed by the postulates of Quantum Mechanics. In Chapter [T] we 
review the mathematical formalism describing quantum measurements and quan- 
tum state transformations. We face problems, such as quantum state preparation 
and repeatability of quantum measurements, which, even though they reach back 
to the beginnings of quantum theory, have been revived and put into a new light 
by the recent developments in the experimental techniques. 

Chapter 121 is concerned with the analysis of quantum channels. Exploiting the 
convex structure of the set of channels, we explicitly single out those that constitute 
the best quantum versions of the intrisically classical copying-, broadcasting-, and 
NOT-gates. We introduce the general theory on which such optimization relies, and 
present group theoretical techniques to analyse the common situation in which 
symmetries of the set of input quantum states, after the action of the channel, 
propagate to the output. This is the framework of covariant channels. It is very 
useful to describe many physical situations and, at the same time, it permits an 
analytical approach. 

Quantum channels are more general to describe changes of quantum states than 
unitary evolutions controlled by Schrodinger's equation. Nonetheless, it is well 
known that every quantum channel is the transformation that a system undergoes 
when unitarily interacting with an auxiliary quantum system — the so-called an- 
cilla — that is discarded after the interaction took place. Actually, this is the only 
way to deterministically realize a non- unitary quantum channel. In Chapter El we 
propose feasible implementations for some of the channels constructed in Chap- 
ter 121 providing the ancillary quantum state and the global unitary interaction. 
This is just a first step towards the experimental realization which remains a far 
more difficult task, however, the setup we propose to optimally copy quantum sys- 
tems, in the case of qubits (i. e. two-levels systems) coincides with the one already 
used in experiments. This is encouraging in view of a possible generalization of 
experimental techniques to higher dimensional quantum systems. 

The thesis ends with Chapter |3] which deals with classical and quantum noise. 
Noise is considered as acting both on the measuring apparata and on the quan- 
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turn states. More specifically, we introduce a (partial) ordering on the convex set 
of measuring devices. It allows us to characterize "clean" devices, namely, those 
which are not affected by quantum and/or classical noise. Interestingly enough, we 
show that such ordering is able to single out von Neumann's observables as "par- 
ticularly nice" measuring apparata. This gives an operational characterization 
for the usually postulated concept of observable. We then focus attention on the 
specific model of noise called decoherence. Decoherence acts destroying quantum 
superpositions, thus making ineffective all quantum improvements on the classical 
approach. On the other side, decoherence possesses also foundational interest since 
it represents the favourite tool to explain the quantum-to-classical transition. The 
process called decoherence is actually a convex set of commuting channels satisfy- 
ing very restrictive properties. Applying techniques described in Chapters El and El 
we provide a method to invert decoherence and restore quantum superpositions by 
a feedback control from the environment. This means that measuring a suitable 
observable of the environment's degrees of freedom, and then performing on the 
system a suitable unitary transformation dependent on the measurement result, it 
is possible to completely cancel the effect of decoherence. 
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Chapter 1 



Quantum Measurements, 
Operations, and Physical Models 

1.1 Classical and quantum events 

Given a finite probabilistic space f2 = {1, . . . , N}, it is possible to define probability 
distributions P = {pi, . . . ,Pn} on Q, where < pi < 1, YliPi — 1- The set of all 
probability distributions on Q, P{fl), is a convex set. It is simple to recognize its 
extremal points as the delta-distributions Pi = 6ij. Such a structure for P(f2) can be 
rephrased saying that P{^1) is a simplex, namely, a convex set whose elements are 
uniquely expressed as a convex combination of extremal points. Random variables 
on are defined as mappings X from into a set of "values" T. Such values 
can be numbers, tensors, or whatever objects. When T is a real vector space, it 
is well-defined the mean value of X : Q — > T, given P e P(il), as X = ^iPiX{i). 
The set of random variables on forms a commutative algebra (under point- wise 
multiplication) . Events are particular random variables where T is the two- values 
set {0, 1}. In the classical case, events form a boolean algebra^: Given two events 
£■1, £'2 : — > {0, 1}, once defined two binary operations A and V as 

E,AE2 = E,-E2, EiW E2 = Ei + E2-Ei-E2, (1.1) 

^ A boolean algebra 58 is a set of elements = {a, 6, c, . . . } satisfying the following properties: 

(i) *B has two binary idcmpotent, commutative, and associative operations, A (logical AND) and 
V (logical or); (ii) *B contains universal boimds and /; (iii) for all a G *B, there exists its 
complementary element a' G 25 such that a Aa' — 0, and a\J a' = I. 
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where "■" is the point-wise multiplication, it is straightforward to verify that all 
properties of a boolean algebra are satisfied. 

Consider now a A^- dimensional complex vector space M' . The analogue of 
probability distributions are N x N density matrices p, i. e. positive semi-definite 
trace-one matrices. The analogue of random variables are N x N hermitian ma- 
trices X. Since random variables, usually called observables, are hermitian, they 
admit a spectral decomposition X = '^jXjUf, where Uf are orthogonal projec- 
tions (of rank greater than one, in case of degeneracy). Density matrices, usually 
called states, define probability distributions over the spectrum of an observable, 
by means of the formula pf{xj) = Tr[pn^]. The mean value of an observable X, 
given a state p, is well-defined as X = Tr[pX]. The non-commutative analogue of 
events are projections Ei = Ef. The set of quantum events E{Jif), called quantum 
logic, has two binary operations A and V defined as 

El A E2 = E1E2, El V E2 = El + E2 — E1E2, (1-2) 

where now the multiplication is the usual (non-commutative) matrix multiplica- 
tion. 

The fundamental differences between the classical model and the quantum 
model are the following (and they are basically equivalent): 

1. the quantum logic is not a boolean algebra, since the distributivity law does 
not hold (beacause of the non-commutativity of the matrix product); 

2. the convex set of states on Jif is not a simplex, but it is strongly convex, 
whence quantum states admit many equivalent ensemble decompositions; 

3. the algebra of observables on Jif is non-commutative. 
See also the introduction paragraphs in p[j and j2j. 

1.2 Notations 

To each quantum system, it is associated a complex separable Hilbert space J^, 
equipped with the inner product {ip\(f)), linear in (p and antilinear in ip, following 
Dirac notation. The set of bounded operators on will be denoted as B{Jif). 
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An operator X is called self-adjoint if it is densely defined and X — on its 
domain^. Self-adjoint operators are called observables and are in correspondence 
with orthogonal resolutions of the identity by means of the formula 

/ + 0O /■ + 0O 

xdU^ix), 7= / dU^{x). (1.3) 
■oo J —oo 

Positive semi-definite trace-one operator p e T+(J^) are called state. We will 
denote the set of states of a system M' as S(Jf ). Since they are all compact opera- 
tors, states can be essentially viewed as infinite density matrices, also in the infinite 
dimensional case, with no relevant differences from the usual finite dimensional set- 
ting. Prom now on, if not otherwise specified, we will deal with finite d-dimensional 
Hilbert spaces isomorphic to C'', for which all linear operators are everywhere de- 
fined, bounded and trace-class, and the self-adjointness coincides with hermiticity. 
Moreover, spectral resolutions are all discrete, \. e. X — Y^- XjHf . 

Composite systems carry a tensor- product Hilbert space. Mi ® M'^ ® ■ ■ - ^ M'n- 
Bounded operators B(J^) form themselves a Hilbert space isomorphic to 

® M' = M'®'^. Once fixed a basis b = for J^, we define the following 

isomorphism between operators in B(J^) and vectors in M"®"^: 

X = ^X,,K)0-| « \X)) ^ ® Ij), (1.4) 

ij ij 

satisfying 

1. {{X\Y)) = Tr[X^Y], i. e. the Hilbert-Schmidt product; 

2. {X (g) Y)\Z)) = \XZY'^)), where Y'^ denotes the transposition with respect 
to the fixed basis b; 

3. Tri[|X))((y|] = X'^Y*, where Y* denotes the complex conjugation with re- 
spect to b; 

4. Tr2[|X))((F|] = Xrt. 

^An operator is called hermitian if its domain is dense in ^ and X C X'^ . In finite dimension 
the two definitions coincide and there is no need to bother with the density of the operator's 
domain. 
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With this notation the state \I /\fd)) is the maximally entangled state on Jf®^: 

-L|,».-L^|,:)«|.>. ,,5) 

Such a state will play a major role in the characterization of quantum devices. 

1.3 Quantum measurements statistics: POVM's 

Given a state p and an observable X = Ylj^j^f^ statistical postulate states 
that: The probabihty of obtaining a result xj within a set A = {xj}jQj is given 
by 



p{xj G A) = Tr 



X 



1.6) 



This means, as we already saw, that a state induces a probability measure pf{xj) = 
TrlpU-j-] over the set of outcomes for a given observable. 

It is clear that, apart from the actual measured value Xj of the observable 
X, the statistics of the outcomes is completely determined by the structure of its 
spectral resolution {nj*-}. In the case of an observable, such H/s are orthogonal 
projections, i. e. HfHf = ^fSy, summing up to the identity, J^j = ^- With a 
little abuse of terminology, from now on we will refer to an observable just as a set 
of orthogonal projections resolving the identity, and to the observable outcomes as 
the indices j's labelling different x/s. 

The concept of observable is generalized by the concept of positive operator- 
valued measure (POVM, for short), which is a set of positive operators P = 
{Pi, P2, . . . , Pn} summing up to the identity J^i^i — ^- Notice that Pj's need 
not to be orthogonal, not even projections, and the number of oucomes of P, i. e. 
its cardinality |P| = N, can be larger than the Hilbert space dimension d. As 
before, also in the case of POVM's, the probability of obtaining the j-th outcome, 
given the system in the state p, is postulated to be p^{j) = Tr[pPj]. 

We call a two-outcomes POVM P = {P, I — P} an effect or, equivalently, a 
property. According to ^ we say that an effect P = {P, I — P} describes a real 
property for the system J^f in the state p, if Tr[pP] = 1. 
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Finally, we introduce here the definition of range^ of a POVM, a concept we 
will extensively use in Chapter HI 

Definition 1.3.1 (POVM range) Given a POVM P = {Pi, P2, . . . , Pn] , its 

range, denoted as Rng(P), is defined to he the convex set of probability distributions 
P = {PijP2, ■ ■ ■ ,Pn} obtained as Pi = Tr[pPj], varying p in all S{J^). 

Remark 1.3.2 Notice that, since p in Definition \1. ,7. il moves around the whole 
quantum states' set, the range of a POVM identifies uniquely the POVM. In other 
words, the correspondence 

P „ Rng(p) (1.7) 

is one-to-one. 

1.4 Quantum operations and instruments 

Since now, we dealt only with the outcomes statistics. However, in order to com- 
pletely describe the measurement statistics we need also to specify the state re- 
duction from prior state p to posterior state pj conditioned by the outcome j. The 
state reduction is nothing but a rule telling us which is the system's state after 
the measurement has been performed and the outcome collected. 

1.4.1 State collapse postulate 

Von Neumann jl] derived the well-known state collapse rule starting from the 
following hypothesis: 

1. the observable to be measured has discrete spectrum and it is non degenerate, 
namely all its eigenspaces are one-dimensional, in formula X = Xi\xi) 

2. the measurement is perfectly repeatable^: Literally from von Neumann's book 
"z/ a physical quantity is measured twice in succession in a system, then we 
get the same value each time" . 

^There is no possibility of confusion between the range of a POVM and the range of an 
operator, being two completely unrelated concepts. 
'^See Section O 
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If such hypotheses are verified, then the system state after the measurement is 

P ' — Pj = (1-8) 

Liiders jH] generalized von Neumann's theorem to degenerate observables, in- 
troducing the postulate of minimum disturbance in the sense that a state, for 
which a property P is real, is left unchanged by a measurement of P. According 
to Liiders' rule, when measuring the observable X = ^ . XiUf the system's state 
after the measurement is 

The interpretation problems to which the state collapse postulate led are be- 
yond the aim of this manuscript. 

1.4.2 Quantum operations 

The appropriate mathematical objects describing a general quantum state change 
are the so-called quantum operations jH]. A quantum operation S, is a completely 
positive trace- non-increasing linear mapping from T'^{Jif) of an input system Jif 
to T+(J^) of an output system Jt^. The map £ : T+(Jf ) T+(J^) is gener- 
ally probabilistic, and the trace Tr[£^(p)] < 1 represents the probability that the 
transformation 

occurs. Deterministic quantum operations, i. e. completely positive trace-preserving 
maps such that Tr[£^(p)] = 1 for all p G S{Jif), are called channels. All quantum 
operations admit the highly non-unique Kraus representation 

£ip) = J2EjpEl (1.11) 

j 

where Ej^s are linear operators from Jif to J^. Nonetheless, it is always possible to 
choose a Kraus representation such that TT[EjEj] = WEiW^Sij] we call it canonical 
Kraus representation. A quantum operation is a channel if and only if its Kraus 
operators satisfy the normalization condition 

Y,E\E, = I. (1.12) 
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Remark 1.4.1 The Liiders' recipe for state change clearly corresponds to a quan- 
tum operation Sj{p) = UfpUf. Notice that the average reduced state p = J2iPi'^)Pi 
can be read as the output of the channel £{p) = UfpUf. 

Every quantum operation S : S{J^) S{Jif) induces naturally a quantum op- 
eration from B(J^) to B(Jf ) by means of the duality relation Tr[£^(p)X] = 
Tr[p£^'^(X)], valid for all p G S{J^) and X G B(^). The map is called the dual 
map, and £'^{Ix) = Ij^f if and only if £^ is a channel. 

Remark 1.4.2 Given a measurement whose outcomes statistics is described by 
means of the POVM P, there exist many different channels associated with P. 
These channels are written as £^{p) = J2i^Fip)' '^'^^^ i^fYi^) = -^i? choos- 
ing between them correspond to assign a particular state reduction rule. 

1.4.3 Instruments 

In the modern formulation of Quantum Mechanics, the most general tool used 
to describe statistical correlations between the outcomes of successive measure- 
ments is given by the notion of (completely positive) instrument, which has been 
introduced by Davies and Lewis [Zj. An instrument is basically a mapping J 
from the set fl of outcomes to the set of quantum operations on S{Jif), such that 
J(lJ^gjj) = '^j^j'3{j) and 3{fl) is a channel. The fundamental result about 
instruments is the following [H] 

Theorem 1.4.3 (Ozawa, 1984) Every statistical measurement theory, consist- 
ing both of outcomes statistics and state reduction rule, can be described by means 
of an appropriate instrument. 

Actually, instruments formalism has been introduced in literature mainly to handle 
the case of continuous outcome space Q, which is described as a standard Borel 
space equipped with a a-algebra 03 (fi). When Q is discrete and subset of M — as 
in our case — technical results become much simpler. For further details on the 
general case see [9j. 

Finally, we define a perfect instrument as an instrument such that J(j) is a pure 
contractive map, i. e. 3{j){p) = MjpM^- , for all j. For example, Liiders instrument 
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in Remark ll.4.1h s a perfect instrument with Mj = 11^. The instrument in Remark 
[1321 is perfect only if £f{p) = ^jp^j. 

1.5 Physical realizations 

Intruments provide both outcomes statistics and state reduction due to a measure- 
ment process. Implicitly, we assume that such a measurement is nondestructive, 
in the sense that the system is left in a state conditioned by the outcome and not, 
for example, absorbed by a counter or a calorimeter. The only reasonable way to 
look for an implementation of a nondestructive measuring process on a quantum 
system is to engineer an indirect measurement scheme. This means that we 
make the system interact with an apparatus and, after some time, we measure 
an observable Y on the apparatus. In formula: 

a(j)(p) = Tr^ {{l^ ® nf) f/(p ® |a)(a|)f/t] , (1.13) 

where {IlJ}j is an orthogonal resolution of Jc/ coming from the diagonalization of 
Y G B(^). Clearly such a procedure gives rise to an instrument, as described in 
the previous Section. Ozawa jH] proved the converse: 

Theorem 1.5.1 (Ozawa, Indirect Measurement, 1984) Every instrument U 
admits an indirect measurement scheme as in Eq. M.l'^) . 

The correspondence is not one-to-one: there are many different — though statisti- 
cally equivalent — indirect measurement schemes producing the same instrument; 
conversely, given the indirect measurement scheme, the resulting instrument is 
unique. 

1.5.1 Levels of description of quantum measurements 

There are basically three ways to describe the statistical aspects of quantum mea- 
surements, depending on the level of details required: 

1. One is interested only in the outcome statistics. Then the maximum gen- 
erality lies in the concept of POVM, as we saw in Section Notice that, 
given the outcome statistics for all quantum states, the POVM is defined 
uniquely — see Remark [1.3.21 
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2. Also the state reduction rule is requested. The notion of instrument encloses 
all possible cases, see Section 11.41 Evidently, many different instruments 
produce the same outcome statistics, i. e. they all correspond to the same 
POVM. 

3. The most detailed description characterizes even the state of the apparatus, 
the physical interaction between the system and the apparatus, and the 
observable to be measured on the apparatus. Clearly the same instrument 
is obtainable by means of different indirect measurement schemes. 

Summarizing, given the outcome statistics, the POVM is uniquely defined. Given 
the POVM, there are many instruments describing it. Similarly, there are many 
indirect measurement schemes realizing a given instrument. The choice between 
different equivalent physical realizations of a measurement process can be made 
according only to "practical" considerations. 

1.5.2 Example: standard coupling 

Consider a discrete observable X = of the system Jif in initial state \ip) 

and let the apparatus system be ^ = L^(]R) in initial state \(f)a)- Now, let and 
interact in such a way that the the observable X couples with the apparatus' 
momentum Pa- This means that the unitary operator is 

U = e-'^^®^\ (1.14) 

The momentum operator Pa is the generator of translations, in the sense that 

e-^^'«^'"0,(x) =0,(x-a;o), (1.15) 

where (j)a{x) = {x\(f>a)- Hence, the initial system+apparatus state \4')^\(pa) evolves 
as 

U\^) ^ <j)a{x) = J2^f\^) ® - fl>^X,), (1.16) 

i 

and, by making assumptions on the value of the coupling constant A and the initial 
state \4>a), if is always possible to obtain functions (paix — HXxi) with (almost) 
disjoint supports. In other words, it is always possible to model an interaction 
between system and apparatus such that the indirect measurement is a position 
measurement on the apparatus — the usual "pointer position" measurement. 
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1.5.3 Example: embedding and optimal phase measure- 
ment 

From the geometrical point of view, the unitary interaction of the system with a 
fixed ancilla state in the indirect measurement scheme simply corresponds 

to a linear (isometrical) embedding of the system ^ into a composite Hilbert 
space M' ® . The measurement on the apparatus then defines a conditional 
expectation from M' ® £^ to M' ^ giving rise to probability and state reduction. 
An embedding into a larger space can always be described by means of an isometry 
V , i. e. a bounded operator such that V'^V = I. If the input system state is p, then 
the embedded state — that is, the system+apparatus state after the interaction — is 
U{p®\a){a\)U^ = VpVl 

In Ref. [inj, we exploited an embedding for single-mode states of the elec- 
tromagnetic field in order to achieve a physical realization of the optimal phase 
measurement. It is well known that the phase of the electromagnetic field does 
not correspond to any self-adjoint operator. Quantum estimation theory [H ITT] 
provides the optimal POVM for the phase measurement in terms of Susskind- 
Glogower operators 

dm = ^\e^'i'){e% l\ili<P) = I, (1.17) 

where |e'"'^) = Yl'^=o^^'^^\'^) ■ The optimal phase measurement outcomes distribu- 
tion is then 

Using the double-ket notation introduced in Section lTT^ consider now the eigen- 
states of the hetherodyne photocurrent Z = a — b'' 

Z\Diz))) = z\Diz))), (1.19) 

where D{z) = e^"^~^*° are the displacement operators, satisfying the completness 
relation 

'^'''D{z))){{D{z)\=I^'. (1.20) 



TT 



The following isometry 



V = ^ [ d'afi\a\)\D{a))){e'^^^% H ^t\fit)\' = - (1-21) 



27r Jc Jo 
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embeds a single-mode state into a two-modes state in such a way that, measuring 
the hetherodyne photocurrent 

p{z) = -Ti[VpV^ \D{z))){{D{z)\] 

\ (1.22) 

= 2 1/(1^1)1' 

one obtains the optimal phase distribution d/i^ as the marginal of ^(2;) on the 
variable = arg z. Notice that here we are performing a joint measurement on 
both modes, not just an indirect measurement on the second mode. However, the 
form of the embedding V provides a natural way to implement the phase POVM 

(HHD. 

1.6 Repeatable measurements 

In Subsection 11.4.11 we introduced the von Neumann-Liiders state collapse prin- 
ciple, derived from the hypothesis of discreteness of spectrum, repeatability, and 
minimum disturbance. In what follows, we derive all the consequences that arise 
from the only hypothesis of repeatability, thus obtaining the most general form of 
a repeatable measurement. See p2| for a detailed derivation. 

First of all, why should we focus on repeatable measurements? Clearly, there 
are a lot of natural measurement schemes which are far from being repeatable, 
think of e. g. a photon counter or a fluorescent screen at the end of a Stern-Gerlach 
apparatus. In the past decades, however, technology of quantum experiments 
improved in such a way that nondestructive measurements on individual atomic 
objects are quite a common task, see e. g. one atom micro-masers and ions traps. 

In the modern formulation of Quantum Mechanics, repeatability hypothesis 
has lost the in-principle relevance it enjoyed in the early foundational books as 
von Neumann's. Nowadays, repeatability is understood just as a property which 
characterizes some particular measurement processes. More precisely, repeatable 
measurements are related to preparation procedures. In fact, preparing a quantum 
system in a particular state means preparing it in a state having some pre-specified 
real property, as defined in Section ESI For example, in order to prepare the pure 
state lip), one may take a collection of quantum systems and perform over them a 
repeatable measurement of the effect described by the POVM {\ip){ip\, I — 
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Of course, the preparation succeeds when the outcome comes out. Then, 

the von Neumann-Liiders state collapse rule tells us that the state of the system 
after measurement is in fact the pure state \ip). In this sense, repeatable mea- 
surements have often been regarded as measurements of observables — projective 
orthogonal resolution of the identity — causing a collapse of the state on one of 
their eigenvectors. 

In ^2] we showed that there exist repeatable measurements which give rise to 
nonorthogonal POVM's and, moreover, which do not even admit any eigenvector, 
that is to say, the reduced state is different at every repetition of the measure- 
ment. This result makes a clear separation between the concepts of repeatability, 
preparation and reality in Quantum Measurement Theory. 

The starting point is the hypothesis of repeatability. A first consequence of 
this is due to Ozawa 0: 

Theorem 1.6.1 (Ozawa, Repeatable Measurements, 1984) An instrument 
satisfies repeatability hypothesis only if it has discrete spectrum. 

Then, perfect^ repeatable instruments are described by a set of contractions {Mj} 



for all j, k and all G J^. The only technical point recalled in the paper is 
that, allowing for infinite dimensional Hilbert spaces, one has also to deal with 
properties of operators such as closeness. In our case, all Mj's are bounded and 
everywhere defined, hence closed [13j. A close operator possesses closed range 
and kernel (the support is always closed since by definition it is the orthogonal 
complement of the kernel) and the Hilbert space M' can hence be decomposed as 



such that Y.i Ml Mi = I and 




(1.23) 



~ Ker(M,) © Supp(M,) ~ Rng(M,) © Rng(Mj)^, 



(1.24) 



The consequences are the following: 



1. All ranges, for different outcomes, must be orthogonal, i. e. 




(1.25) 



^For the definition of perfect instruments, see Subsection 11.4.31 
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2. All ranges must be contained in respective supports, i. e. 

Rng(Mi) C Supp(Mi), Mi. (1.26) 

3. All Mj's satisfy the condition 

M/M,|R,g(M,)=/Rng(Af,). (1.27) 

Now, the fundamental difference between operators on finite and infinite Hilbert 
spaces is that, in finite dimension, Supp(X) ~ Rng(X) always, while, in the infi- 
nite dimensional case, one can have Supp(X) C Rng(X) or, viceversa, Rng(X) C 
Supp(X), strictly. This holds basically because in the infinite dimensional case 
there exist proper subspaces with the same dimension as the whole Hilbert space 
M^. This observation lead us to the following: 

Theorem 1.6.2 For finite dimensional systems, only observables admit repeat- 
able measurement schemes, and the system state collapses according to the von 
N eumann-Luders rule M.f^) . 

So the finite dimensional case describes precisely what one usually expects about 
the structure of repeatable measurements. It is nonetheless possible to construct 
a simple example in infinite dimension, enclosing all counter-intuitive features of 
the infinite dimensional case. Let us consider a two-outcomes POVM: 



Po=p|0)(0| + X^|2j + l)(2j + l|, 

oo (1-28) 
Pi = (l-p)|0)(0|+^|2j + 2)(2j + 2|. 

j=0 

Notice that P = {Pq,P{\ is a nonorthogonal measurement. We can describe such 
a POVM by means of the following instrument: 

oo 

Mo = v^|l)(0| + |2(J + 1) + l)(2j + 1|, mIm, = Po, 

(1-29) 

Ml = VT^|2)(0| + J2 |2(J + 1) + 2)(2j + 2|, M|Mi = Pi, 

3=0 
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in the sense that, got the i-th outcome, the state changes as MipMj . Repeata- 
bility hypothesis can be simply checked. Analysing the structure of scheme fll.29j) 
one recognizes a unilateral-shift behaviour of the kind 5*1^) = |n + 1). Actually, 
this unilateral-shift structure is a general feature of nonorthogonal repeatable mea- 
surements. Since S does not admit any eigenvector, analogously the scheme ()1.29|) 
changes the system state at every repetition of the measurement and there are 
no states which are left untouched by such a scheme. In other words, in infinite 
dimensional systems there exist repeatable measurements which cannot satisfy 
minimum disturbance hypothesis, even in principle, and hence cannot be viewed 
as preparation procedures. 



Chapter 2 

Characterization and 
Optimization of Quantum Devices 

In order to handle information encoded on quantum states we need to engineer 
astonishingly precise and accurate devices since the least loss of control in manip- 
ulating quantum systems can lead to extremely detrimental effects on the whole 
process. The theoretical investigation is the starting ground in designing such op- 
timal quantum devices. This Chapter is devoted, first, to giving a complete and 
tractable characterization of quantum channels, second, to exploiting such char- 
acterization to single out optimal devices according to particular figures of merit 
that we will introduce and explain from time to time. 

The basic assumption we will adopt is to consider input quantum states be- 
longing to sets obeying some symmetry constraints — i. e. satisfying invariance 
properties under the action of some groups of transformations. Moreover, we will 
choose figures of merit conforming in a natural way to the same symmetry con- 
straints. These two conditions lead to the very well established mathematical 
framework of covariant channels, for which the characterization simplifies, mak- 
ing explicit calculations analytically solvable. Actually, covariant channels form 
convex sets whose structure is (in some cases) known and optimal devices lie on 
the border of such sets. In this way, the problem resorts to a semi-definite finear 
program. 

In particular, we will focus on channels optimally approximating the impos- 
sible tasks of copying, broadcasting, and performing NOT on unknown quantum 
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states. The symmetries we wiU deal with are universal symmetry (invariance un- 
der the action of SU((i)), phase-rotations symmetry (invariance under the action 
of U(l)^'^), and invariance under the group of permutations. 



2.1 Choi-Jamiolkowski isomorphism 

A useful tool to characterize quantum channels in finite dimensional systems is 
the Choi-Jamiolkowski |^ |S] isomorphism — one-to-one correspondence — 
between channels E : S{J^) S(J^) and positive operators Re on J^^Jf defined 
as follows: 

Re = {S0I)\I)){{I\ ^8 = Tr^ [(/® p^) Re] , (2.1) 

where T is the identity map on S{J^), \I)) = \i) ® \i) is the maximally en- 
tangled (non normalized) vector in ,Jif ® Jif, and denotes the transposition 
with respect to the fixed basis used to write |/)). Different Kraus representa- 
tions for S{p) = ^.iEipE\ correspond to different ensemble representations for 
Re = X]j \Ei)){{Ei\, the canonical^ being the diagonalizing one. Trace-preservation 
constraint J^iM-^i — rewrites as Tr^ \Re\ = ^.X'- 

Choi-Jamiolkowski isomorphism ()2.1|1 turns out to be very useful in describing 
covariant channels. In the following Section we shall recall some basic notions 
about group theory. 



2.2 Group-theoretical techniques 
2.2.1 Elements of group theory 

A unitary (projective) representation on M' of the group G is a homomorphism 
G 9 (? 1-^ f/(, G B(J^), with \Jg unitary operator, such that the composition law is 
preserved: 

UgUH = iu{g,h)Ugh- (2.2) 



^See Subsection HXl 
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The cocycle u}{g, h) is a phase, i. e. \uj{g, h)\ = 1, for ah g,h E G, and it satisfies 
the relations 

uj{gh, k)u{g, h) = uj{g, hk)uj{h, k) 

A unitary representation is called irreducible (UIR) if there are no proper sub- 
spaces of M' left invariant by the action of all its elements. Two irreducible 
representations and C/^ of G on and J^, respectively, are called equivalent 
if there exists a unitary T : J^i ^ ^ such that TU] = U'^T, for all g e G. The 
fundamental result concerning UIR's of a group is the following: 

Lemma 2.2.1 (Schur) Let and U'^ be two UIR of G on M\ and M2, respec- 
tively. Let B : Mx — > M2 a (bounded) operator such that: 

BUl = Up, (2.4) 

for all g & G. Then: 

1. U^ and U"^ equivalent =^ B oc T; 

2. U^ and U^ inequivalent =^ B — 0. 

Remark 2.2.2 (Abelian groups) From Schur Lemma simply follows that fact 
that, if the group G is abelian, nemely gig2 = g2gi for all gi,g2 G G, then all 
its UIR's are one- dimensional. In fact, all Ug's must be proportional to the same 
unitary operator T and they are all simultaneously diagonalizable, hence reducible 
on direct sums of one- dimensional invariant subspaces. 

2.2.2 Invariant operators and covariant channels 

Let Wg a reducible unitary representation of G on J^f. Then J^f can be decom- 
posed into a direct sum of minimal invariant subspaces: 

^~0^. (2.5) 

i 

Each -Jifi supports one UIR of G. Some UIR's can be equivalent or inequivalent. 
Let us group equivalent UIR's under an index /x labelling different equivalence 
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classes, and let an additional index span UIR's among the same /z-th equivalence 
class. Since equivalent UIR are supported by isomorphic subspaces, i. e. M'i^ ~ 
~ M'p, for all iyi^ip, in the same /x-th class, we can rewrite the decomposition 
(Q as 

Jr~0Jr^®C'^^ (2.6) 

where c?^ is the cardinality (degeneracy) of the /i-th equivalence class. Decompo- 
sition ()2.(ij) is usually called Wedderburn's decomposition T7], the spaces are 
called representation spaces, and the spaces multiplicity spaces. Then, the 
following decomposition for the representation Wg holds 

W9 = ®Wg^^h^. (2.7) 

With Eq. ()2.7|1 at hand, it is simple to derive the form of an operator B, invariant 
under the action of the reducible representation Wg, i. e. 

W^gBWg = B, 'ige G. (2.8) 

Since the above implies [B, Wg] = 0, for all g, then: 

5 = 0F®5,^^, (2.9) 

where Bd^ is an operator on C^** . In other words, the operator i? is in a block-form 
since it cannot connect inequivalent representations and can act non-trivially only 
on multiplicity spaces of the representation Wg. This is precisely what is contained 
in the Schur's Lemma [2.2. II 

Now, consider a family of quantum states F C S{J(f) that is invariant^ under 
the action of a group G, namely UgpU^ G F for all G G and all p G F. The 
group, and then the family F, can be discrete as well as continuous. A channel 
S : F ^ S(eJ^) is said to be covariant under the action of the group G if 

SiUgpUD = Vg£{p)Vl yg G G, (2.10) 

where Ug and Vg are two generally reducible unitary representations of G on 

and respectively. In a sense, the channel S is "transparent" with repect to the 

^Notice that this requirement is weaker than requiring that F is the orbit of a single seed state 
under the action of G. 
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action of the group G and the image of the invariant family F is another invariant 
family S{F). Using Choi-Jamiolkowski isomorphism ()2.H) . the above covariance 
condition for S rewrites as an invariance condition for Rs namely, 

[Re,Vg®U;] = 0, WgeG, (2.11) 

where, as usual, the complex conjugate is with respect to the basis used to write 
Rs. Decomposing ^ (g) = 0^ (g) one gets: 

Re = ^I''0Rd,, (2.12) 

with positive blocks Rd^. 

Another direct consequence of Eq. ()2.9|1 is the form of a group- averaged oper- 
ator, namely 

{X)g= [ dgU^XUl [ dg = l. (2.13) 
Jg Jg 

Clearly, (X)g is invariant, whence, if the representation Ug of G decomposes the 

Hilbert space as Jf = 0^ C^^*, it can be written as 

Notice that Tr^[X] is a short-hand notation for Tr^ [P^XP^], where P^ is the 
projection of Jif onto Jf^ C^^'. 

2.2.3 Example: §I[J((i)-covariance 

A typical §I[J(c/)-covariance, also known as universal covariance, for short U- 
covariance, is that under the representation of many input and output copies, 
namely when Jif = {C^)^^ and = (C'^)®^^ with Ug = W^^ and Vg = W®^ . 
Here, Wg is the defining representation of SU((i), and invariance condition ()2.11|1 
reads: 

[Rs, PVf ^ ® {W;f^] = 0. (2.15) 

The general Wedderburn's decomposition for such a representation is very com- 
plicated and channels satisfying covariance ()2.15|1 will be studied with a somewhat 
different approach, see Subsection 12.4.11 Nonetheless, there are two situations in 
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which universal covariance can be conveniently faced using Rs machinery. The 
first situation is when Ug = Wg and Vg = W*. This is the case in which we are 
requiring a controvariance condition: 

S{WgpW^^) = W;£{p)Wl (2.16) 

The invariance condition reads \Rs^ IV®^] = 0, which implies Rs = rsP^^^ +rAP^\ 

(2) (2) 

where Pg and P\ are respectively the projections onto the totally symmetric 
and the totally antisymmetric subspaces of Jif'^'^. We will analyze this case in 
Subsection 

The second situation is when d = 2, namely when we deal with qubits. First 
of all, in this case the two representations Wg and W* are equivalent, since 
W* = OyWgOy Heuce the Wedderburn's decomposition for Wf^^ ® (W;)®^ 
is the same as for wf'^^^^^^ which the well-known Clebsch-Gordan series ^H] for 
the defining representation of §0(2)^: 

Af/2 Af/2 

"^©^ ' " ' (2.17) 

M/2 7V/2 j+l 

-00 ^J^C'^^^C'^S 
i=jo i=«o J=\j-i\ 

where jo; are equal to or 1/2 if M, are even or odd, respectively, and 

2 j + 1 / M \ 

'^^■-M/2+j+iU/2-jr ^^-^^^ 

We will analyze this case in Subsection 12.4.31 
2.2.4 Example: U(l)-covariance 

The defining representation of U(l) is simply a phase e**^ G C. In higher dimen- 
sions, we can impose either 'phase- covariance jT], that is, 

U^ = e''^^, N = n\n)(n\, n = 0,...,d-l, (2.19) 



•^Rigorously speaking, this is not the Wedderburn's decomposition since different can 
support equivalent representations. See Subsection l2.4.3l 
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useful to model systems driven by a Hamiltonian with equally spaced energy 
levels, as the harmonic oscillator Hamiltonian, or multi-phase covariance, that 
is, covariance under a unitary representation of the d-fold direct product group 
U(l) X ■■■ X U(l): 

d-l 

U^ = Y^e''l^-\n){n\, 0„G[O,27r[, (2.20) 

n=0 

where (p = {</>„} is a vector of d independent phases. Notice that one of such phases 
is actually an overall phase and can be disregarded: for a (i- dimensional system 
we then have {d — 1) effective phase-degree of freedom. In the following we shall 
adopt multi-phase covariance, and, where there is no possibility of confusion, we 
shall interchange the terms phase- covariance and multi-phase covariance. Notice 
that, in the case of qubits, the two concepts coincide. 

Also phase-covariance is typically applied to many copies of input and output 
(say and M, respectively). For qubits the representation U^'^ decomposes as 
(see Eq. flTTjl l 

N/2 

f;«iV^0,.0ji%j^^^ (2.21) 
l=lo 

where = Yln=-iMh''^){hM is the angular momentum component along rota- 
tion axis, say 2;-axis, in the / representation. As in the universal case — U(l) is a 
subgroup of SU(2), actually — dealing with two-dimensional systems allows us to 
handle the complete Wedderburn's decomposition and work in full generality, even 
with mixed sates (see Subsection I2.5.3|l . 

In higher dimensional systems, we shall restrict ourselves to pure input states. 
This implies that the many-copies input state lives actually in the totally symmet- 
ric subspace"^ = (C^)! ^. Moreover, optimal map will be found to have output 
supported in = (C'^)f*^. Now, a convenient way to decompose the composite 
space ® in the Wedderburn's form 0^ (g) C^** , is the following: 

JT ® ~ ® jr, (2.22) 

{m-i} 

'''This is true only for many-copies pure input states . Indeed, a many-copies mixed input 
state p®^ is generally non symmetric. 
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where {rrii} is a multi-index such that YliJ^i = M — N. Invariant subspaces are 
clearly one-dimensional, since the group is abelian, and equivalence classes are 
spanned by^: 

=^{mO®^ = Span(|{mi + ni})® |K})| . (2.23) 

I. J {Ui} 

In the above equation, {rii} is a multi-index such that Yli^i ~ The vectors 
llrij}) are defined as: 



IK}) = 




(2.24) 



where 11:^^ are permutations of the systems. In other words, are totally 

symmetric normalized states, whose occupation numbers are denoted by the multi- 
index {rii}. Clearly, by varying {rii} over all possible values < rij < A^, the 
set spans all input space J^. Analogous arguments hold for the vectors 

\{mi + rii}) in That the decomposition using |{mj + rii}) ® is useful 

to identify the block structure of a multi-phase covariant channel is clear noticing 
that 

jj^M ^ (f/*)«A^|{^^ + m}) ® IK}) = e^E.-.0.|{^. + n,}) ® IK}), (2.25) 

for all possible choice of {rii}. We'll make use of this decomposition in Subsections 
IZOandEEl 

2.2.5 Example: permutation group invariance 

Most channels of physical interest act on input states which are indeed "many- 
identical-copies states". This is the case, for example, of estimation channels, 
which optimally reconstruct an unknown input state by performing measurements 
on A^ copies of it. Analogously, when the task is distributing quantum information 
to M users, typically one requires that the reduced state is the same for each user. 

^We consider here only maximally degenerate equivalence classes, namely, equivalence classes 
whose degeneracy equals the dimension of the input Hilbert space (C*)®^. For example, the 
vector (g) |0)®^ supports an irrep but it cannot be written as in Eq. (|2.23() . In Subsection 

12.5.11 we will see how this constraint indeed does not cause a loss of generality. 
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Both situations can be described by saying that input and/or output states are 
actually permutation invariant states. In formula: 

£{p) = ^(nf pnf ) = nf ^(p)nf , Vr, a, (2.26) 

where U!^ and 11^^ are (real^) representations of the input and output spaces 
permutations, respectively. When both properties are satisfied, the operator 
must equivalently satisfies the following invariance condition: 

[i?^,nf ®nf] = 0. (2.27) 

Notice that such an invariance property is stronger than that in Eq. ()2.11|) since 
it implies both conditions 

[Rs, nf ® J®^] = [Re, I^^' ® nf ] = 0, (2.28) 

for all a, r, hence in particular for a = t. The fundamental tool that comes now 
at hand is the so-called Schur-Weyl duality between permutation group represen- 
tations on qubits and the defining representation of SU(2). The duality relation 
tells that decomposes (C^)*^*^ precisely as W^'^ , namely, 

M/2 

(C2)®M _ (^2i+l ^ ^d, ^ ^2.29) 

j=jo 

but with exchanged role for the spaces. Explicitly, C"^^^^ is now the multiplicity 
space and C^^ the representation space. In turns, from Eq. ()2.9|1 . Schur-Weyl 
duality gives the form of a generic permutation invariant operator X on (C^)®^: 

M/2 

[X, nf ] = ^ X = ® J,^, , (2.30) 

j=jo 

where is an operator on C^-'^^. 
Decomposition of many-copies qubit states 

As an application of Schur-Weyl duality, let's consider the decomposition of the 
many-copies qubit states p®^. This decomposition has been first given in Ref. pUj . 
For a complete and detailed proof see Ref. PT] . 
^Representations of permutations are always real. 
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Indeed, such many-copies states are invariant under permutations of single 
qubit systems. The state p®^ admits then a decomposition as in Eq. ()2.3()jl . 
explicitly 

N/2 N/2 I , . „ 



^ ©E (2.31) 

l=lo n=-l ^ ^ 



where, as usual, p = (/ + rfc ■ cr)/'i, ||fc|| = 1 and \l.,n) are the eigenvectors of 
the angular momentum along fe, namely jf'. Notice that Eq. ()2.3H1 exhibits a 
singularity for r = 1 due to the particular rearrangement of terms. However, the 
limit for r ^ 1 exists finite, as it can be seen from the equivalent expression 

. 2\ N/2-1 I I X i+n 



©^ E 137 \l,nWM®W (2.32) 



2.3 Optimization in a covariant setting 

Let us consider a family F = {po} of quantum states of the input system In 
most cases of physical interest, such a family is invariant under the action of a 
unitary representation Ug on of a group G, in formula: 

UgPeUl = pgio) e F, Vp G F, \/g e G. (2.33) 

On such a family of states we are concerned about a particular mapping 7V1 of F 
onto another family F' = {ag} of states of the output system invariant under 
the action of another unitary representation Vg of the same group G. The mapping 
Ai can be completely general, even physically non allowable. Let Ai be covariant, 
namely, M{pe) = erg. 

Whatever Ai is, we introduce a physical channel S and a merit function ^, 
depending on 6 and Ai, such that ^[S{pg), ag] = ^{6) achieves its maximum when 
S{pg) = cxg. In other words, ^ quantifies how well the channel S approximates the 
mapping At. Assuming transitive action of G on F, that is, 

ye, 3geG: e = g{eo) for a fixed ^o, (2.34) 

a further natural requirement is the invariance property of ^: 

mdo)) = ^ [SiUgPeoUD, Vgag,V^^] = ^ [V^^SiUgpg.UDVg, ag,] = m)- (2-35) 
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Then the function to be maximized is the average score 



^= / 5^ [^(Ps{eo)).^<;(eo)] dg = d{Oo)- 
Jg 



(2.36) 



The basic point is that, if the optimum average score is reached by some channel 
S, it is always possible to achieve the optimum also by a covariant channel S, 
namely such that S{UgpUl) = VgS{p)Vgl Indeed, from Eqs. and it 

turns out that'' 



where we defined 



^(Pgido 



(2.37) 



V9^iP9ieo))Vg dg. 



(2.38) 



It is simple to verify that [R^, Vg ® U*] = 0, namely, S is covariant and, by 
construction, it achieves the optimal average score ^. 

Hence, in the following we can restrict the optimization procedure to covariant 
channels, which form a convex set. By introducing appropriate convex merit func- 
tions, we can moreover search for the optimum channel within the border of the 
convex set, since convex functions defined on convex sets achieve their extremal 
values on the border. In the cases in which we are able to characterize extremal 
covariant channels, we can then explicitly single out channels optimizing the given 
merit function. 



2.4 Universally covariant channels 

Universal covariance means, in literature, covariance under the action of the group 
SU((i). Invariant families of states contain states with fixed spectrum: the most 
usual choice is to restrict the analysis to the set of pure states. Given a channel 
S : S{J^) S(J^), universal covariance reads S{UgpU^) = VgS{p)V^ where Ug 

^We also assume ^ linear in the 1. h. s. slot. 
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and Vg are unitary representations of §U(c?) on and respectively. In the case 
of pure input states {ip), we will consider as merit function the fidelity, namely, 

m\^)m, i0)(0i] = Tr s{\^)m] = tt [{\<j>){<j>\ ® \^){^n rs] ■, (2.39) 

in the case of mixed input (qubit) states p = (/ + raz)/2, we will consider the 
purity (the Bloch vector length^), namely, 

d[Sip), z] = Tr [a, S{p)] = Tr [(a, ® p*) Re] . (2.40) 

It is clear from the form of score functions ()2.39|1 and ()2.4()|1 that both are convex 
(linear) in Rs and invariant (see Eq. ()2.35|l ). 

2.4.1 Optimal universal cloning 

In this Subsection we shall basically review Ref. j22] using Choi-Jamiolkowski iso- 
morphism. We can't thoroughly apply the formalism we developed in the previous 
Sections because a closed form for Wedderburn's decomposition of f/^^ represen- 
tation of SU((i) is very complicated. We will follow a somewhat alternative path, 
finding a particular map maximizing the score function and satisfying covariance 
and trace-preservation conditions^. 

Quantum cloning of an unknown state p®^ —>■ p®*^, M > A^, is impossible 
|23j . Much literature has then been devoted to searching for optimal physical 
approximations of impossible ideal cloning [24 . Two basics assumptions are made 
in order to make calculations treatable: such optimal machines should work equally 
well on all input states, and input states should be pure p = ip = The 
natural framework to work within is then the universal covariance. The score 
function is taken to be the fidelity between the actual output of the approximation 
map Ci^ilj®^) and the ideal output ip®^^ . In terms of the Rc operator: 

J[C(V'^'^), ^^'"] = Tr [(^«*^ ® {rf") Rc] , (2.41) 

where Rq, in order to satisfy universal covariance of C, is such that 

[Rc, U'f' ® (t/;)®^] = 0. (2.42) 

^Sometimes the purity is defined to be proportional to the square of the Bloch vector length: 
Tr[p2] ^ (l + r2)/2. 

^For uniqueness proof see Ref. 
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Let P. 



(iV) 



P. 



(TV) 



be the projection over the totally symmetric subspace of 



the input system M' = (C^)^^ . Since ip^^Pf '' = ^j^^ , we have that 



< Tr 



P< 



(TV) 



The channel C is universally covariant, whence, from Eq. ()2.14j) . 



C (Pf )) = I <lgUf^C (Pf (f/t)«^ 



Tr 



CP, 



iP^^ + O, 



(2.43) 



(2.44) 



where O collects all other contributions coming from partially symmetric/antisymmetric 
invariant subspaces, and d[M] = is the dimension of the totally symmetric 

subspace. Actually, terms in O does not contribute to the fidelity since ip^'^ is a 
symmetric state^°, hence, w. 1. o. g., we write 



CP, 



,(JV) 



d[N] (M) 

and obtain the following upper bound for the score function: 

d[N] 



d[M] 



One can easily verify that the positive operator 



Rc. 



d[N] 
d[M] 



P, 



(A/) 



®TVv 



r(g)N 



y(M) 



(2.45) 



(2.46) 



(2.47) 



is invariant, properly normalized to trace-preservation^^, and saturates the bound 
()2.4(jj) . With a little abuse of notation, we denoted with |J®^)) the non-normalized 
maximally entangled vector in (C^)®^^ 

d-l 

|/®^)) = |i7v)®Ki)®---® Ktv), (2.48) 

■ -n' ' ' ^ ' 

ii,...,ijv— u (^d-j^iv (cd)(8JV 

such that 

Tr [(^^^ ® (^*)^^) |/^^))((/«^|] = Tt[^Y = 1, (2.49) 

since i/) is pure. From Choi-Jamiolkowski inverse formula, one can verify that the 
action of the optimal universal cloning is as given in Ref. [^, that is, 

(2.50) 



^"Here it is crucial that p = i/; is pure. Otherwise p"^^ could also have non-null components 
on partially symmetrized/antisymmetrized subspaces. 



11 



In the sense that Tr j(f[Rc] = IjPs 
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2.4.2 Optimal universal NOT-gate 

Another unphysical mapping with a naturally emerging covariant structure is the 
NOT-gate. In this Subsection we shall derive, following Ref. [23], the optimal phys- 
ical approximation of the ideal quantum-NOT. Actually, in Subsection 13.2. ![ we 
shall also show that optimal cloning and optimal NOT are intimately related. 

Let us consider a d- dimensional system M' described by the pure state = 
When ci = 2, it makes sense to consider the NOT-gate, which, generalizing 
the classical mapping — > 1 and 1^0, sends an unknown pure state to its 
unique orthogonal complement. Such orthogonal complement, a part from a fixed 
unitary transformation, is the transposition of the input state. This fact explains 
why perfect NOT-gate is not physical, since transposition is the simplest example 
of positive transformation that is not completely positive. In [2E] the case d = 2 
is addressed and the optimal universal approximation is worked out. Here we 
generalize the result for all finite dimensions and pure input states. 

First of all, it is clear that for d > 2 the orthogonal complement of a pure state 
is not uniquely defined. Hence we shall construct the map T approximating the 
transposition, which, on the contrary, is uniquely defined — once fixed a basis in 
Jif. Universal covariance for a channel whose output transforms as the transposed 
input, that is, T{UgpU^) = U*T{p)Uj , reads, as usual, as an invariance property 
for Rr: 

[Rr,U;^U;] = 0. (2.51) 

The unitary representation (f/*)®^ of SU((i) decomposes the space J^®^ into the 
irreducible totally symmetric and totally antisymmetric subspaces, J^g''^ and Jf^^ 
respectively. Hence Rt = rsP^ + taPa K where P^^ : J^^^ Jf^l are orthog- 
onal projections. 

The covariant score function is taken to be the fidelity Ti[^*T {ip)], as always 
when dealing with pure states. From the form of Rq-'- 

^ = Ti[{rT^Rr]=rs, (2.52) 

and rs has to be maximized consistently with trace-preservation condition 
Tr^fr^P^^^] = /. Noticing that pf^ = (J®^ + S)/2, where S is the swap-operator 

(2) 

between the two spaces, its partial trace is easily computed as TY,^^[rsPs ] = 
Irs{d + l)/2. The optimal universal approximation of the transposition map is 
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then uniquely described by 

Rr = ^PP, (2.53) 

and it achieves optimal fidelity 

d = 2/{d+l). (2.54) 

Remarkably, such value for the fidelity equals the fidelity of optimal state es- 
timation over one copy [2^. This means that, even if the optimization has been 
performed in a general setting, the resulting channel T, that is optimal and unique, 
corresponds to nothing but a trivial measure- and-prepare scheme. In other words, 
optimal universal transposition can simply be achieved by performing the optimal 
state estimation over one copy — the input copy — and then preparing the trans- 
posed of the estimated state. This aspect is usually referred to as classicality 
of the channel. We will see in Subsection 12.5.21 that, in the case of multi-phase 
covariant transposition, this classical limit can be breached. 

2.4.3 Universal qubit superbroadcasting 

Broadcasting of quantum states is a generalization of cloning, in the sense that 
given an unknown input state p G S(^), the broadcasting machine B is allowed to 
return a generally entangled output S G S^Jf^"^) such that Tri[S] = Tr2[S] = p. 
In j2H] it's been proved that it is not possible to broadcast with the same channel 
two noncommuting quantum states. This result is generally referred to as the no- 
broadcasting theorem. Actually, the proof holds only for single-copy input state; 
allowing for multiple-copies input, it is possible to construct a channel broadcast- 
ing a whole invariant family of states. Moreover, considering as merit function the 
Bloch vector length (in the case of qubits, see Eq. (jTlOJ), the optimal broadcasting 
channel actually purifies the input state, in the sense that the single-site reduced 
output commutes with the input (hence their Bloch vectors are parallel) being at 
the same time purer (with longer Bloch vector) than the input. We will refer to 
such a broadcasting-purifying gate as the superbroadcaster |,29]. Of course, the su- 
perbroadcaster can be made a "perfect" broadcaster by appropriately mixing the 
output state with the maximally chaotic state 1/2 (this procedure simply corre- 
sponds to a depolarizing channel isotropically shrinking the Bloch vector towards 
the center of the Bloch sphere). 
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In what follows we will explicitly derive such an optimal superbroadcasting 
machine by thoroughly using group-theoretical techniques of Section I2I21 

Permutation invariance and universal covariance 

We consider a map B taking copies of an unknown qubit state p to a global 
output state of M > qubits. A first natural requirement is that each final 
user receives the same reduced output state^^. This fact, along with the obvious 
permutation invariance of the input p®^, leads to a Choi-Jamiolkowski operator 
that must satifsy the following invariance property (see Eq. ()2.27|l ): 

[nf ®nf,/?e] =0, Va,r, (2.55) 

where 11^ and Ilf^ are (real) representations of the permutation group of the M 
output and the A^ input systems, respectively. From Eq. ()2.30|1 the form of Rb 
follows 

A//2 N/2 

= ® ® Id, , (2.56) 

3=30 l=lo 

where Rji is an operator on C^-'^^ ® C^'^^ and dj and di are the Clebsch-Gordan 
multiplicities given in Eq. ()2.18j) . 

Eq. ()2.56|) takes into account only permutation invariance of input and out- 
put states: it can then be further specialized to different group-covariances. In 
this Subsection we will deal with SI[J(2) covariance (see Subsection 12.5.31 for U(l) 
covariance). According to Subsection I2.2.3[ since W* = ayWgay, such covariance 



condition rewrites as 



a 



0, where Wg is the defining representation 



of §U(2) and St3 = (/®^^ ® a^^)Rts{I®^^ ® a®^). Hence splits as 

M/2 N/2 j+l 

^s =00 (2.57) 

j=jo l=lo J=\j-l\ 

where P^'i is the orthogonal projection of the space C^-''^^ C^''^^ onto the J- 

representation and satisfies the simple properties: 

2 7-1- 1 2 7 -1-1 

TT[Pf,] = 2J + 1, TtAP/,] = Tt,[P/,] = (2-58) 



^^This requirement alone could not give rise to permutation invariant output states. However, 
it is possible to prove that one can always find an optimal map satisfying this property, see 



Ref. EH. 
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Classification of extremal points 



Since Sb has to be positive, all s^^ are positive real numbers and trace-preservation 



condition Tr.x[S'e] = reads 

N/2 M/2 j+l 

21 + 1 



Tr^[^^] = © E E ^^m^I^^-^^ ® = I''"'- (2-59) 



1=10 3=30 J=\j-l\ 

The latter is equivalent to 

27+1 

E E »^Ai^ = l. V/. (2.60) 

3=30 J=\j-l\ 

To single out optimal maps, here we adopt the Bloch vector length merit function 
()2.4()j) . This is a linear merit function, thus optimal maps lie on the border of the 
convex set of covariant channels described by Sq operators. The problem is how 
to characterize extremal Sq operators compatible with complete positivity and 
trace-preservation constraints. Since Ss is diagonal in indeces j and J, extremal 
Sb operators are classified by functions j = ji and J = Ji, leading to the following 
expression for extremal 5*^ operators 

l=lo 

Optimization 

We now feed the input state p®^ into the channel. Since we are working in a 
universally covariant setting, we can write, w. 1. o. g., p = (/ + raz)/2, that is, an 
input state with Bloch vector along z-axis. The global output state S is 

S = Tr^[J«^^ ® (^,pV,)«^ Sb] = Tr^^[/«*^ ® p«^ Sb] , (2.62) 

where p denotes the NOT of p, corresponding to the inversion r —r (or, equiva- 
lently, r± r^). By means of the decomposition ()2.31|1 for p®^, we get 



l=lo n=-l 



(2.63) 
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p (r) 




0.2 0.4 0.6 0.8 1 

Figure 2.1: The plot shows the behaviour of the scahng factor p^'^+i(r) versus r, 
for ranges from 10 to 100 in steps of 10, in the universal case. Notice that there 
is a wide range of values of r such that p^'^+i(r) > 1. 

From the form of the map, it is guaranteed that S is permutation invariant. Hence 
it makes sense to speak about the reduced state a = Ttm-i[^], regardless of which 
particular reduced state. In |2I] it is shown that [o"^,^] = 0, namely, the reduced 
state Bloch vector is along 2;-axis. The merit function is then 

d = TT[{a,®l''-') S] =r', (2.64) 

where r' is the Bloch vector length of a. After a lengthy calculation (see ^^), 
the optimal channel turns out to be the one with ji = M/2 and Ji = M/2 — I, 
regardless of the number of input copies and of the spectrum of p. The optimal 
superbroadcasting achieves the following scaling factor p^'^^(r) = r'/r: 

^ ' l=lo n=~l 
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The two hmiting cases are M = N + 1 



P 



,N,N+l 



(r) 



{N+l)r 



N + 3 



^ ^ l=lo n=-l ^ 



1 — r 



1 + r 



) 



n 



(2.66) 



and M = oo 




r\ 4 / ^ l + l ^ \l + r 

^ ^ l=lo n=-l ^ 




n 



(2.67) 



By plotting scaling factors for different values of and M, it turns out that, 
in the universal case, superbroadcasting first emerges for = 4 (in Subsection 
I2.5.3l we will see that, in the phase-covariant case, superbroadcasting first emerges 
for = 3). Quite surprisingly, for a sufficiently large number of input copies 
[N > 6) it is possible to superbroadcast quantum states even to an infinte number 
of receivers. In Fig. 12. II there are the plots oi p^'^~^^{r) for 10 < < 100 in steps 
of 10. Notice that for r 1 all curves go below one: indeed optimal universal 
cloning of pure states never achieves fidelity one, see Eq. ()2.46|) . 

A compact way to describe the performances of the optimal superbroadcaster 
is to introduce the parameter r*, implicitly defined by the equation 



Clearly, r* actually depends on A^ and M. By the monotonicity of p{r), for r < 
r* there is superbroadcasting. Hence, r* > means that superbroadcasting is 
possible. As we already said, for A^ > 6, r* > for all M. For A^ = 5, r* > 
for M < 21. For A^ = 4, r* > for M < 7. Moreover, as A^ and M get 
closer, r* — 1, as expected. In Fig. 12.21 there are the plots of 1 — r*{N,M), for 
M = A^ + 1 and M = oo. With good approximation, the two curves have power 
law 1 - r*(A^, A^ + 1) oc 2/A^2 ^ _ ^) ^ 

2.5 Phase-covariant channels 

Multi-phase rotations in d dimensions, see Eq. ()2.20|1 . obviously form normal sub- 
groups of SU((i). In other words, multi-phase covariance group is "smaller" than 
SlLJ((i) and, consequently, multi-phase invariant families of states contain "less 




(2.68) 
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Figure 2.2: Logarithmic plot of (1 — r*) versus in the universal case. The upper 
line refers to the case M = oo and shows a behaviour like The lower line is 

for M = iV + 1 and goes like 2/N'^. 
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states" than universally invariant families. Actually, multi-phase invariant fami- 
lies directly generalize in higher dimension the idea of the equator of the qubits 
Bloch sphere^^. Quite intuitively then, optimization in a multi-phase covariant 
setting should generally achieve better performances than the analogous univer- 
sal optimization, since the group is smaller and leaves margin to sharperly tune 
free parameters. In what follows we will consider the same examples of the previ- 
ous Section (cloning, NOT-gate, and superbroadcasting) in a multi-phase covariant 
framework and we will compare the results. 

2.5.1 Optimal phase-covariant cloning 

The task is to optimally approximate the impossible cloning transformation ip®^ — > 
^ where ip is an unknown pure state belonging to a family invariant under the 
transitive action of the multi-phase group, whose defining representation is 

d-l 

f/^= |0)(0| + ^e^'^"|n)(n| (2.69) 

n=l 

(with respect to Eq. ()2.20jl here we put 00 = 0, since an overall phase is irrelevant). 
As before, since we deal with pure states, the input space is considered to be 
the totally symmetric subspace (C'^)f^. Analogously, the output space is J(f = 
(C'^)f'^. Invariant figures of merit are the usual (global) fidelity 

[c{i^r).i^r] = Tr [cii^r) , (2.70) 

and the single-site fidelity 

[TrM^i [C{i^r )\ , ^0] = Tr [C(^o'^^) (^o ® , (2.71) 

where ■j/'q = dr^/'^ YltZl K) is a fixed state whose orbit spans all possible input states 
family. We will adopt 5^^, nonetheless, in Ref. j^I] we proved that multi- phase 
covariant cloning maps optimizing single-site fidelity optimize global fidelity as 
well. Clearly, it is understood that the channel C satisfies the covariance property 

[Rc, ® (f/;)n = 0' (2-72) 

"'^^This idea can be made more rigorous noticing that, when d + 1 mutually unbiased basis 
can be written, d of them are connected by multi-phase rotations, as it happens for qubits. See 

Ref. isni- 
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so that Eq. ()2.70|1 makes sense. Last condition leads to the following form for Rc 

Rc = Yl E 4<}!k'}Iw.} + k})(w.} + k}i®im^ (2-73) 

{m,}K},K'} 

where we used the compact notation defined in Eq. ()2.24|) . As usual, Rc has 
to be positive, in order to guarantee complete positivity of C, and satisfy trace- 
preservation condition Tr j^[_Rc] = 

After lengthy calculations (see Ref. |31j); the optimal multi-phase covariant 
cloning machine is found to be the one described by the positive rank-one operator 

Rc= \{^^ + k})^{< + k}\®\{n^}){{<}l (2.74) 

K},K} 

where A; is a positive integer such that + ~ ^i hence equal to {M — N)/d. 

Optimal single-site fidelity is 



A^! {ni + k + l){nj + k + l) 



d Md^+^ ^ ^ n^\...ni\...n.\...\\ (n^ + l)(n,- + 1) 



(2.75) 



which for A^ = 1 simplifies to 



1 (rf-l)(M + rf-l) 

In Fig. 12. HI there are the plots versus M of optimal 1 — > M single-site fidelity in 
the cases of multi-phase covariant and universal cloning for c? = 5. Multi-phase 
covariant cloning achieves better fidelity than the universal one, as expected. 

Notice that our analysis is not completely general because of the restricting 
relation that must hold between input and output number of quantum systems 
involved 

M = N + kd, keN. (2.77) 

However it is the most general result on multi-phase covariant cloning machines 
described in the literature by now. 



2.5.2 Optimal phase-covariant NOT-gate 

The multi-phase covariant approach to approximate the NOT-gate is one of the 
examples in which the performances improvement, with respect to the universal 
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case, is more apparent. The transformation we consider is the NOT-gate tjj ^ ijj* 
for pure ci-dimensional states belonging to a multi-phase invariant family spanned 
as before by the multi-phase rotations group applied to a fixed seed state 
tpo = d"^/"^ 'Ylii V)- The covariant figure of merit is the fidelity 

i?[r(^o),^o]=Tr[r(^o)^o], (2.78) 

since ipQ = ipo (with the appropriate choice of basis) . The channel T must satisfy 
the covariance property 

[Rr,U;^U;]=0. (2.79) 

The group is abelian so that all irreps are one-dimensional. Equivalence classes 
with respective characters are classified in Table 12.11 The Rr operator then splits 
into a direct-sum 

i? = 0i?n0% (2.80) 

i i>j 

of 1 X 1 blocks Rii acting on Span{|z) ® and 2x2 blocks Rij acting on 
Span{|^)®|j),|j)®|i)}. 
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|0)® |1),|1)®|0) 


g-«01 







Table 2.1: Equivalence classes and respective characters of irreducible one- 
dimensional representations of (f/^)®^- 

In Ref |S21 there is the complete derivation of the final form of optimal Rr 
operator as 

i>j 

where bij > are matrix elements of a null-diagonal symmetric bistochastic^^ 
matrix. For d = 2,3 this constraint suffices to single out a unique optimal T, since 
the only null-diagonal symmetric bistochastic matrix for c? = 2 is 




(2.82) 



and for d = 3 





( ° 


1/2 


l/2\ 


{b^A = 


1/2 





1/2 




\l/2 


1/2 


/ 



Already for d = 4, there are two free parameters < j9i < 1 and < p2 < 1 — Pi 



A matrix is called bistochastic if all its rows and columns entries sum up to one |33|. 
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in defining a null-diagonal symmetric bistochastic matrix 



( 

V\ 

V2 
Vl - V\2 



Pi P2 1 - Pl2\ 

1 - Pl2 P2 

1 - Pl2 pi 



P2 



The achieved optimal fidelity is 



Pi 



Pl2 = P1+P2- 



(2.84) 



/ 



(2.85) 



strictly greater than in the universal case ()2.54|) . for all d. Moreover, it is interesting 
to notice that 2/dis also greater than the fidelity of optimal multi-phase estimation 
over one copy, derived in Ref. to be {2d — l)/d'^. This means that, contrarily to 
the universal case for which the optimal NOT-gate is classical (see final remarks in 
Subsection I2.4.2|) . the multi-phase covariant analogue breaches the classical limit. 
The result is particularly striking in the case of qubits for which it is not possible to 
perfectly estimate the phase with finite resources, while it is possible to perfectly — 
with unit fidelity — transpose an unknown pure equatorial state by means of a fixed 
unitary transformation. 



2.5.3 Phase-covariant qubit superbroadcasting 

In the phase-covariant version of superbroadcasting, we specialize the permutation 
invariant form ()2.56p imposing the further constraint 

[RB,U^^^{U;r'']=0. (2.86) 

Let us now suppose that input states lie on an equator of the Bloch sphere, say 
xy-plane. Then, Urp are precisely rotations along ^-axis, namely 

U^ = e'^"% (2.87) 

and invariance condition ()2.86|) rewrites as 

= 0, Vj,/, (2.88) 



•J. tU) -j. r(i 
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where jT* = Y^n=-i'^\^i^)^^^'^\ angular momentum component along z- 

axis in the / representation. A convenient way to write operators Rji satisfying 
Eq. (jT88jl is the following: 

I I j-i 

^^•'= E E E r^^'Jj,n + k){j,n' + k\(^\l,n){l,n'\, (2.89) 

n=—l n'=—l k=l—j 

when j > I, and 

j j l~3 

^^•'=E E Y.'''^,rn',k\3^^)^3^^'\®V^^ + k){l.m' + kl (2.90) 

m=—j m'=~j k=j—l 

when j < I, both expressions exhibiting similar structure as in Eq. ()2.73|) . Notice 
that there are two more running indeces with respect to the universal case ()2.57|1 . 
While the index n' in Eq. ()2.89|1 simply allows for off-diagonal contributions, 
the index k labelling equivalence classes is related to the direction of the reduced 
output state Bloch vector, as we will see. In particular we will show that, in order 
to get an equatorial output, operators Rji have to be symmetric in k, in the sense 
that ri^y^j^ = 

Classification of extremal points and /c-symmetry 

Trace-preservation now reads 

EE<n,.^^ = l' (2-91) 

j k 

and, analogously to the universal case, the fact that r^'^ ^ > and Rji operators 
are diagonal with respect to indices j's and fc's implies that extremal points are 
classified by functions 

3=3u k = ki. (2.92) 

Equivalently, extremal Rji are proportional to correlation matrices^^ since they are 
positive matrices with diagonal entries r^''^ all equal to l/dj^ (see Eq. ()2.9H) ). and 
extremal correlation matrices are known in literature j^. In particular, rank-one 
correlation matrices are extremal, hence rank-one operators Rji are extremal. 

^^Correlation matrices are positive semi-definite matrices with diagonal entries all equal to 
one. 
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Figure 2.4: Schematic sketch of the fc-symmetrization procedure. 



In order to further simplify the general form of Rb in Eqs. fl2.89|l and ()2.9()j] . 
we now impose on the single-site reduced output state the following additional 
constraint 



Tr 



M-l 



I 

2' 



(2.93) 



We will see that constraint ()2.93|) . on one hand, does not cause a loss of generality 
since it does not affect optimality, and, on the other, clarifies the geometrical 
interpretation we mentioned about fc-indexed degrees of freedom of phase-covariant 
broadcasting maps. In fact we have (for explicit calculation see Ref. 











TrM-1 


K2OJ 


= TrA,f_i 


Tr.^ 



Re 



Since J^iC^^ + l)di = 2^, the only condition for Eq. is that 



^2.94) 



0. 



(2.95) 



In a sense, index k labels a "tilt" of the reduced output state Bloch vector with 
respect to the equatorial plane. Our requirement is then a "null tilt-requirement" , 
or, in other words, an "equatorial output state-requirement", and it can always 
be achieved by equally mixing two extremal maps — generally losing extremality — , 
the first labelled by a function k = ki, the second by k = —ki. We will refer to 
such a property as fc-symmetry property of bradcasting maps and we showed that 
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/c-symmetry property is equivalent to the property of mapping equatorial states 
to equatorial states. Notice that a A;-symmetric map is such that r^'^, ^ = r^'^, 
The strategy to obtain broadcasting maps optimizing the reduced output state 
Bloch vector length is then to search for optimal maps within extremal maps and, 
once found the best one, to force /c-symmetry on it. The procedure is shown in 
Fig. 12.41 On the left there are the equatorial mixed input state p and the single- 
site reduced output p'^. Suppose such an output comes from an extremal map 
Bk described by r^'^, ^ elements. Notice that, by covariance, the projection of 
p'^ onto the equator is parallel with p. Consider now another map B^ki whose 
elements are equal to f^'„, ^ = r^'„, _f^. Clearly, B-k is a proper channel obeying 
all covariance and extremality constraints as Bk- The output of B-k is in sketched 
in the middle figure as p'_^. In order to have an equatorial output, we mix Bk and 
B^k obtaining B = {Bk + i3_fe)/2 whose output p' = (p'^ + p'_^)/2, by linearity, lies 
on the equator (see the picture on the right). Of course B is no more extremal, 
by construction. However, the figure of merit we are considering, namely, the 
length of the projection of the output Bloch vector onto the original one, does not 
change. In this sense, imposing fc-symmetry does not affect optimality. Moreover, 
it is possible to prove that the fc-symmetrized output p' has higher fidelity with 
the input p (see Ref. |2I])) than the tilted p'^ and p'_i^. 

Optimization 

In Ref |2I] it is proved that the channel optimizing the merit function 

d = TT[{a^^l''-') E], (2.96) 

for x-oriented input states p = (/ + rax)/2, has ji = M/2 for all /, and ki = 0, 
for M — even, while ki = ±1/2 for M — N odd. Hence, for M — N even the 
optimal superbroadcaster is already /c-symmetrized, whereas for M — N odd we 
must equally mix the channels coming from ki = 1/2 and ki = —1/2. In both cases, 
^nn' ki ~ l/'^ip ^'y ^- 6'^'^' structure of the map B depends only 

on the parity of M — A^, and not on the spectrum of p. 

For M — N even, the optimum scaling factor p^'^^{r) is given by 

N/2 N/2 I 



\ / i^i^ n=-l L ^ / -I n,n+l 

(2.97) 



n,n+l ' 

n,n+l 
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p (r) 




Figure 2.5: The plot shows the behaviour of the scahng factor p^'^+i(r) versus r, 
for ranges from 4 to 100 in steps of 8, in the phase-covariant case. Notice that 
there is a wide range of values of r such that p^'^+^(^r) > 1. 
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Figure 2.6: Logarithmic plot of (1 — r*) versus in the phase-covariant case. The 
upper line refers to the case M = oo and shows a behaviour like 1/2N. The lower 
line is for M = iV + 1 and goes like 2/{3N^). 



while, for M — N odd, it is 

^ ^ l=lo n=-l L \ / 

(2.98) 

In Fig. 12.51 there are the plots of p^'^+i(r) for 4 < < 100 in steps of 8. 
As in the universal case, all curves, for r — > 1, go below one: indeed optimal 
phase-covariant cloning of pure states never achieves fidelity one, see Eq. ()2.75|) . 
However, it is possible to see that phase-covariant superbroadcasting is more effi- 
cient than the universal one: superbroadcasting first emerges for A^ = 3 (A^ = 4 in 
the universal case) and achieves larger values of p^^^{r) for all A^, M, and r. 

In Fig. ESI there are the plots of 1 - r*(A^, M), for M = A^ + 1 and M = oo, as 
done for the universal superbroadcaster. With good approximation, the two curves 
have power law 2/3N'^ and 1/2N, respectively, namely they go to zero faster than 
in the universal C3iSG, clS expected. 



r^..n.-\-^ 



n+l/2,n+3/2 



Chapter 3 

Realization of Quantum Devices 



In the previous Chapter we exphcitly wrote quantum operations coming out from 
an optimization procedure in a covariant setting. We gave such channels in terms 
of their Choi-Jamiolkowski operators 1^2.1^ . However, Choi-Jamiolkowski repre- 
sentation for quantum channels, even if very useful in dealing with semi-definite 
programming problems, turns out to be quite far from giving the physical "recipe" 
needed to realize the channel in a laboratory. In the following we will describe how 
to unitarily implement a given quantum channel, in terms of a unitary interaction 
between the system and an ancilla. In the first Section, we will provide, following 
Refs. |36| 1^ . a general method to work out a physical setting realizing a given 
channel. In the second part of the Chapter, we will show how this procedure works 
in the case of some of the channels discussed in Chapter |2l 



3.1 Unitary dilations of a channel 

Let us given a channel S : S{Jif) S{M'Y . The task of this Section is to find an 
ancilla system an ancilla pure state |0), and a unitary operator U on ® 
such that 

^(p) = Tr^[[/(p®|0)(0|)f/t], (3.1) 

for all p e S{je). 

^Here, for simplicity we disregard the case of channels from states on a system to states 
on another system J^, e. g. the cloning channel from S{J^'^^) to 5{J^^^). This case can be 
taken into account, see Ref. [37] for a more general approach. 
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3.1.1 Stinespring dilation 

Given a channel the Stinespring representation ^H] (V, of £^ is a kind of 
"purification" of the channel i. e. 

£{p) = Tr^ [VpV^] , (3.2) 

where V is an isometry, i. e. V'^V = I, from Jif to ® The Stine- 
spring representation is usually given for the dual channel (see Subsection ll.4.2j] 
: B(^) ^ B(^) as 

S-{0) = V^O(^I^)V. (3.3) 

Let £{p) = J2i^iPM be a Kraus representation for S. Consider now the 
operators from to ® ^ defined as Ei Cg> where |i) belongs to a set of 
orthonormal vectors in The only trivial condition ^ must satisfy is dim^ > 
^{Ei}. Then, the sum 

V = J2E^®\^) (3.4) 

i 

is an isometry, since V^V = Y^i-^l-^i ~ ^.-^^ ^^"^ realizes the channel 8 as in 
Eq. jH. 

Remark 3.1.1 Notice that we did not make any assumption on the particular 
choice for the Kraus representation {Ei\ used to construct the Stinespring isometry 
V in Eq. When {E^} is the canonical Kraus decomposition and d'ms^ = 

^{Ei}, we will refer to such V as the canonical Stinespring representation for S, 
which clearly is the one minimizing the ancillary resources, i. e. the dimension of 
the ancilla system, needed to physically implement the channel. 

3.1.2 Unitary dilation 

From Stinespring form ()3.2j) the existence of a unitary interaction U between Jif 
and ^ realizing the channel S is apparent, since every isometry V from to 
M' ® £^ can obviously be written as 

V = \J{I^®\^)), (3.5) 

where ?7 is a suitable unitary operator on ® and |0) is a fixed normalized 
state of . Now, |0) is precisely the ancilla state such that 

£:(p) = Tr^[f/(p®|0)(0|)f/t]. (3.6) 
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While the existence of a realization U for every channel is a well-established 
fact in the literature [HI EH], the problem of giving explicitly such interaction for 
a given channel can be very difficult. The general procedure given in Ref. [HZI 
basically relies on a repeated Gram-Schmidt orthonormalizing algorithm applied 
to the column vectors of the Stinespring isometry V. In this way we are able to 
find additional dimJf x (dim^/ — 1) orthonormal vectors to append to V, com- 
pleting it to a square matrix whose column vectors form an orthonormal basis for 
the composite system ® i. e. to a unitary operator on ® £/. In the 
following Section, we will see that, in some fortunate cases, the channels optimized 
in Chapter |2l admit very simple Stinespring isometrics, allowing us to explicitly 
write unitary operators realizing such channels in dimension d. 

3.2 Explicit realizations 

3.2.1 Universal NOT and cloning gates 

Let us start from the optimal universal NOT-gate T derived in Subsection l2.4.2[ The 
channel is completely described by the positive operator Rr in Eq. ()2.53|1 . In order 
to write T in its Stinespring form, we first have to obtain a Kraus decomposition 
for T. This can be done by expanding Rr (see Section 12. ip as 



d-l 




d+1 



1 




1 



d-l 



(3.7) 



d-l 



J2 \ML))m 



rS I 

mn I ? 



where 




1 



{\m){n\ + \n){m\). 
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One possible Kraus decomposition is then given by 

r(^) = J2 MLi^M^n- (3.9) 

m,n=0 

A Stinespring isometry V such that T{ip) = Tr^ [y-?/'!^"''] is then^ 

d-l 

V= M^^®\mn)U, (3.10) 

m,n=0 

where we chose = as ancilla system. Summarizing, we wrote the optimal 
NOT-gate T by means of an isometry V embedding the input system Jif into a 
composite tripartite system Jif ® Jif (S) Jif, in which the last two spaces represent 
the ancilla. 

Tracing VipV^ over the last two spaces, we get the channel T. What happens 
if we trace over different combinations of spaces? In fact, all three spaces are the 
same and there is no reason to consider one of them as the preferred system and 
the remaining ones as ancillae. Actually, tracing VipV'^ over the first space, one 
obtains 

Tn [VW^] = ^pf\l®^)pf\ (3.11) 

namely, the optimal 1 — 2 universal cloning for pure states (see Eq. I2.50p . This 
means that universal 1-^2 cloning and universal NOT-gate are intimately related 
and contextually appear on different branches (spaces) of the same physical setting. 
Such a coincidence has been experimentally exploited for qubits in Ref. [IHl and 
theoretically analyzed and interpreted in generic dimension in Ref. |3T] . Moreover, 
it is possible to prove that Trs [^t/'I^"'"] optimally approximate the transformation 
if) ^ %p* ® %p for pure states. Notice that the cloning map is basis independent, 
whilst the transposition map depends on the choice of the basis, which is reflected 
by the choice of the particular Stinespring isometry V . 

In Ref. jinj it is possible to find the explicit calculation deriving a unitary 
interaction and an ancilla state realizing at the same time optimal approximations 

^Notice that this Stinespring isometry is not the one minimizing ancillary resources. In fact, 
it comes from a Kraus decomposition which is not the canonical one, since the orthogonality 
condition, Tr[MfjM^^ = for all {ij} ^ {kl}, does not hold. However, as we will see in the 
following, this realization allows a very intriguing physical interpretation, see Ref. |41|. 
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of universal cloning and transposition. The unitary operator U on (C ) is 



{p\{q\ - {q\{p\ 



(3.12) 



p=0 



p,q=0 
p<q 



V2 - ^ V2 

where the three sets of isometrics {V^,p}, l^^,?"*!, and from Jif to Jif' 



are defined as 



d-1 



k=0 



k®p)\k®p){k®p\, 

^ d-1 

= ^^\k){\k®p)\k®q) + \k®q)\k®p)){k®q\. 

^ d-1 

= - \k®q)\k®p)){k®q\ 



(3.13) 



Preparing the ancilla state as 



i-l 



d + 



(3.14) 



r=0 



the following identity holds 

U{ij ® \4>)){{4>\)U^ = V^V\ (3.15) 

namely, the operator U in Eq. ()3.12p together with the ancilla state |(/))) in Eq. ()3.14|) 
provide a unitary dilation of the Stinespring isometry V in Eq. ()3.10|) . realizing op- 
timal universal 1 — > 2 cloning as well as optimal universal transposition, depending 
on what we trace out after the interaction. 

In the case (i = 2, we obtain the network model for universal qubit cloning of 
Ref. with 

/l 000000 o\ 
1 
1 
1 
1 
1 
1 
\0 0001000/ 
and |0)) = ^(2|O)|O) + |O)|1) + |1)|O)). 



(3.16) 
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3.2.2 Phase-covariant cloning and economical maps 

In Subsection 12.5.11 we obtained the channel optimally achieving the multi-phase 
covariant N ^ M cloning transformation. The optimal channel has been de- 
scribed, as usual, by giving the corresponding Rq operator in Eq. ()2.74j) . In the 
analyzed cases, i. e. when M = N + kd, where G N and d is the dimension of 
the single copy system, Rc enjoys the relevent property of being rank-one. This 
implies that its canonical Kraus representation contains only one operator, and, to 
satisfy trace-preservation constraint ()1.12|) . such an operator has to be an isome- 
try. Therefore, the optimal multi-phase covariant N ^ M cloning machine Cn,m, 
for M = N + kd, admits the very simple expression 

Cv.mI^^"^) = Vr''V\ (3.17) 

where V : Jif^^ — > j^®*^ is an isometry acting as follows 

V\M) = \{n, + k}), (3.18) 

using the compact notation introduced in Eq. ()2.24p . 

This kind of isometric optimal channels attracted attention in the recent lit- 
erature as economical transformations [ISll^E^, in the sense that, in order to 
physically implement them, there is no need of discarding additional resources, i. e. 
ancillae. In fact, from the point of view of Stinespring representation, multi-phase 
covariant cloning is realizable as 

C^^^ (^«^) = U (^^^ ® |0)(0r(*^-^)) U\ (3.19) 

namely, with respect to Eq. ()3.6|) . there is no partial trace, and the resources needed 
are just the (M — A^) blank copies where convariantly distribute the information 
contained in ip^'^ . 

3.2.3 Phase-covariant NOT-gate 

The optimal multi-phase conjugation map has been derived in Subsection l2.5.2l to 
be 

Rr = J2^^Jm + mmi + m, (3.20) 

i>j 
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where bij > are matrix elements of a null-diagonal symmetric bistochastic (NSB) 
matrix. In this case the map, for > 2, is not unitary or isometric as in the case 
of phase-covariant cloning. Moreover, the fact that in dimension d > 4 there 
exists a whole set of equally optimal maps — in one-to-one correspondence with 
NSB matrices — makes the problem of finding a physical realization much more 
difficult than in the two examples treated before, where the optimal map was 
unique. There are basically two paths one can follow: the first is to search for 
the most efficient realization, i. e. the one minimizing ancillary resources (in 
this case we will tipically single out one particular optimal phase-conjugation map 
achievable using less resources than the others) ; the second is to search for the most 
flexible realization, i. e. the one that spans as many as possible optimal maps by 
appropriately varying the "program" ancilla state and/or is more robust against 
noise (this second kind of realization will clearly require a higher dimensional 
ancilla system to encode a "fault-tolerant" program). 

A good point to start with is the study of the structure of the set of optimal 
phase-conjugation channel, or, equivalently, of the set of NSB matrices. Such 
matrices form a convex set^. On the other hand, every bistochastic matrix is a 
convex combination of permutation matrices — this is the content of the Birkhoff 
theorem |3S|- The null-diagonal and symmetry constraints, however, force the 
convex set of NSB matrices to be strictly contained into the convex polyhedron of 
bistochastic matrices. This fact causes the extremal NSB matrices to eventually 
lie strictly inside the set of bistochastic matrices, generally preventing them from 
being permutations. 

The geometrical study of the set of NSB matrices and its extremal points 
can shed some light on the unusual feature that there exist different "equally 
optimal" maps. The problem arises for dimension at least d = 4. In this case the 
decomposition of the matrix {bij} in Eq. ()2.84|) into extremal components is 



+ P2 



/O 1 0\ 

10 

1 

\0 1 0/ 



/O 1 0\ 

1 

10 

\0 1 0/ 



+ P3 



/O 1\ 

10 

10 

Vi 0/ 



(3.21) 



^This is because their raws and columns are probability distributions 
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where pi,P2,P3 > and pi + P2 + Ps = 1- A natural question is now which optimal 
maps can be achieved with minimal resources. 

More explicitly, for c? = 4, we define three unitaries Ui, U2 and f/3 on C'^ ® 

as 

(^10 T32 \ I T20 T31 \ / T30 T21 \ 

, f/2 = ' ^3 = , (3.22) 

T32 Tio I \T3i T20 I \T2i T30 / 

where Tij = \ + Each of them realizes an extremal optimal multi-phase 

conjugation map (corresponding to = 1 in Eq. ()S.2H) for a given k), namely 

^(P) = E Bi!;%P^^, = Trjf/fc (p ® |0)(0U) Ul], (3.23) 
i>j 

where |0)(0|a is a fixed qubit ancilla state. Hence extremal phase- conjugation maps 
for d = 4 can be achieved with just a control qubit. Notice that the ancilla must 
not necessarily be in a pure state, and the optimal map is equivalently achieved 
for diagonal mixed ancilla state a|0)(0|a + /?|l)(l|a. By adding a control qutrit, 
we can now choose among any of the optimal maps using the controUed-unitary 
operator on ® ® C'^ 

U = Ui® |0)(0| + U2® |1)(1| + f/3 ® |2)(2|. (3.24) 

Any optimal multi-phase conjugation map can now be written as 

%{p) = TTa,b [U (p ® |0)(0U ® ab) U^] (3.25) 

where cxb is a generic density matrix on C^. By superimposing or mixing the three 
orthogonal states {|0), |2)} of the qutrit we control the weights pi,p2,p3 in Eq. 
()3.2H1 via the diagonal entries of ab- In other words, using a 6-dimensional ancilla 
it is possible to span the whole set of optimal maps. 
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Eqs. fl3.22|l - ()3.25|l can be generalized for higher even dimensions^, with 




= 1, . . . , (i — 1, 



i,j=0 



d-l 



U = J2Uk0\k){k\, 



(3.26) 



k=l 



Tr(p) =Tr4f/fc(p®|0)(OU)f/2], 
T,(p) = Tr,,, [U (p® |0)(OU®or,) f/t] 



where UkS are unitary operators acting on C"'/^, t/ is a control-unitary op- 
erator on C''/^ ® C^~^, |0)(0|a is a fixed ((i/2)-dimensional pure state, and 
is a generic {d — l)-dimensional density matrix. The minimum dimension of the 
ancilla space required to unitarily realize an optimal phase covariant transposition 
map is d/2, generalizing the result for d = 4, for which just a qubit is needed (see 
Eq. ()3.23|) ). On the other hand, to span the whole optimal maps set one needs a 
{d — l)d/ 2- dimensional ancilla. 

Finally, notice that realization of multi-phase covariant transposition generally 
needs much less resources than realization of universal transposition: the minimum 
dimension d/2 of the ancilla space in the phase covariant case has to be compared 
with the dimension d'^ required in the universal case ()3.14|) . 



^The case of odd dimensions is much more complicated and will not be analysed here. The 
problem with odd dimensions is that extremal points of the convex set of NSB matrices are not 
permutations. Hence Birkhoff theorem cannot be applied. 
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Chapter 4 



The Role of Noise in Quantum 
Processes 



In the previous Chapters we saw how to optimize transformations over quantum 
systems and how to reahze them by means of physical interactions. Of course, 
processing of quantum systems requires a very high level of control during all 
steps of the experiment. On the other hand, noise — in the sense of uncontrollable 
interactions of the system with the sorroundings — is not always and completely 
avoidable: the only thing the experimenter can do is to reduce it in order to reach 
the desired level of confidence. This can be done by trying to directly control the 
environment, e. g. forcing it into a cavity, or by engineering states, interactions 
and measuring apparata robust with respect to the adopted model of noise. 

In the following Sections we will deal with noise on measuring apparata and on 
states. While in the first part (review of Ref. [IH]) we will face very general models 
of noise — basically, all non-unitary completely positive maps — in the second part 
(review of Ref. |17j) we will focus on decoherence of quantum states, proposing 
a novel correcting scheme retrieving classical information that the decoherence 
process made leak into the environment and exploiting such information to undo 
the noise. 
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data processing 









(A) 



(B) 



Figure 4.1: There are two ways of processing POVM's: (A) the postprocessing of 
the output data and (B) the preprocessing of the input states by a quantum chan- 
nel. The postprocessing is purely classical, whilst the preprocessing is quantum. 



4.1 Clean POVM's 



Let us given a general apparatus performing a measurement on quantum states. 
We know that the most general way to mathematically model it is by means of a 
POVM P, see Section [1.31 Let us now think for a while we don't know how the 
apparatus P works. It could be noisy at the input gate, that is, quantum states 
undergo some uncontrolled transformation £ before being measured, and/or noisy 
at the output, the outcomes being, let's say, shuffled before being read by the 
experimenter. The two situations are depicted in Fig. 14.11 The question is the 
following: Do we have any condition on P that allows us to a priori afflrm that P 
is "clean", i. e. without noisy processing at the input and/or the output? 

Clearly the point can be viewed from the complementary point of view: What 
kind of processings are possible on a given POVM? How does the apparatus change 
after such processings are performed? 



4.1.1 Postprocessing of output data 

The most general postprocessing of a POVM outcomes is a shuffling of with con- 
ditional probability p{i\i) > 0, corresponding to the mapping 



Q: = Y.pmP3^ 

3 



(4.1) 
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where = 1, Vj. To visualize the shuffling (jUTJ, it is useful to think to 

the POVM as a column of operators and p{i\j) as a column-stochastic matrix^ 



fQi\ 

Q2 

\QnJ 



/j9(l|l) p{l\2) ■■■ p{l\m)\ 



p(2|l) p(2|2) 



p{2\m) 



\p{n\l) p{n\2) ■■■ p{n\m) J 



(4.2) 



\PmJ 



Notice that postprocessing generally does not require that P and Q have the same 
cardinality. Relevant examples of postprocessing are: 

1. identification of two outcomes, e. g. j and k are identified with the same 
outcome /, corresponding to p{n\j) = p{n\k) = 6in] 

2. permutation vr of outcomes, corresponding to p{'n'{j)\k) = Sjk- 

When two POVM's P and Q are connected by a mapping of the form ()4.1|) 
for some conditional probability p{i\j) we will write P >~p Q, and say that the 
POVM P is cleaner under postprocessing — postprocessing cleaner, for short — than 
the POVM Q. It is possible to prove that the relation >~p is a pseudo-ordering, 
hence an equivalence relation under postprocessing can be defined as follows 



Definition 4.1.1 The POVM's P and Q are postprocessing equivalent — in sym- 
bols P Q — if and only if both relations P Q and Q >~pP hold. 

We are now in position to define cleanness under postprocessing, namely 

Definition 4.1.2 A POVMP is postprocessing clean if and only if for any POVM 
Q such that Q >~p P, then also P >~p Q holds, namely P Q. 

The complete characterization of cleanness under postprocessing (classical) is given 
by the following theorem (see Refs. liSj) 

Theorem 4.1.3 (postprocessing) A POVM P is postprocessing clean if and 
only if it is rank-one. 

-^That is, a matrix of positive numbers such that all its columns' entries sum up to one. 
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This means that if a POVM is rank-one, we are sure that it does not have a hidden 
noisy postprocessing at the output. Viceversa, the Theorem says that it is not 
possible to obtain the statistics of a rank-one POVM by classically postprocessing 
the outcomes of a higher rank POVM. 

4.1.2 Preprocessing of input states 

A preprocessing £ of input states induces naturally a dual channel E"^ acting on 
the POVM itself, as seen in Subsection 11.4.21 Hence, we will write 

P ^ Q (4.3) 

and say that the POVM P is preprocessing cleaner than Q, if and only if there 
exists a channel — i. e. a trace-preserving, completely positive map £ — such that 
Qi = £'^{Pi), Wi, or, equivalently Q = £''(P), for short. It is possible to prove that 
the relation :^ is a pseudo-ordering, hence an equivalence relation under prepro- 
cessing can be defined as follows 

Definition 4.1.4 The POVM's P and Q are preprocessing equivalent — in sym- 
bols P ^ Q — if and only if both relations P Q and Q P hold. 

We are now in position to define cleanness under preprocessing, namely 

Definition 4.1.5 A POVMP is preprocessing clean if and only if for any POVM 
Q such that Q :^ P, then also P Q holds, namely P ~ Q. 

From the above definition, it turns out that a POVM is preprocessing clean 
if and only if, whenever a noisy preprocessing acts, its action on the POVM can 
be perfectly inverted. Now, a result by Wigner tells that a channel admits an 
inverse channel (i. e. it is physically invertible^) if and only if such a channel is 
actually unitary. A question arises: Does cleanness property define an interesting 
structure in the set of POVM's? Or will we find that invertible preprocessings are 
just unitary (i. e. trivial) preprocessings? Generally, this is not the case, because 

^There exist channels that are invertible in the sense that they define a one-to-one correspon- 
dence between states, but their inverse mappings are not channels. This is the case, for example, 
of the isotropic depolarizing channel p i-^ pp + {1 — p)I/d. In Ref. ^Bj we actually derived, as a 
corollary, that one-to-one channels either are unitary or their inverse map is not even positive. 
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we want the action of the noise to be invertible only on a fixed POVM, not on all 
B(^). 

For qubits, however, preprocessing-equivalence coincides with unitary-equivalence 

Theorem 4.1.6 (qubits) For two-level systems P ~ Q if and only if there exists 
a unitary operator U such that P = [/"''Qf/. 

In higher dimensions the counterexample is given implicitly by the following 
Theorem regarding effects (two-outcomes POVM's, see Section IT!^ 

Theorem 4.1.7 (effects) Let P = {P,/-P} and Q = {Q,I-Q] he two effects. 
Then P ~ Q z/ and only if Xm{P) = Xm{Q) o^nd Xm{P) = ^m{Q), where Xm{0) 
(\m{0)) is the minimum (maximum) eigenvalue ofO. 

Since necessary and sufficient condition for preprocessing-equivalence of two effects 
is that they have the same spectral width, ragardless of the spectrum itself, it 
is clear that there exist preprocessing-equivalent effects which are not unitarily 
equivalent (otherwise they should have the same spectrum as a whole). It is also 
clear that, for dimension c? = 2, the spectrum is completely determined by the 
spectral width, whence unitary-equivalence. 

Besides effects, the other case in which we have a complete characterization of 
preprocessing clean POVM's is the following 

Theorem 4.1.8 (observables) For number of outcomes n < d, the set of pre- 
processing clean POVM's coincides with the set of observables. 

This result is interesting since it provides an operational approach, alternative to 
the axiomatic one given by von Neumann, to define what are the observables in 
quantum theory. Here, just by introducing the cleanness pseudo-ordering, we sin- 
gled out the set of observables, as the only clean POVM's with numer of outcomes 
less or equal to the dimension of the Hilbert space — in this sense, they are the only 
clean "classical" POVM's. 

When the numbers of outcomes gets larger than the dimension of the Hilbert 
space, the structure introduced by the preprocessing pseudo-ordering on the convex 
set of POVM's becomes more complicated, and we have just partial results. For 
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example, we can prove that rank-one POVM's are not only postprocessing clean, 
but also preprocessing clean 

Theorem 4.1.9 (rank-one) Rank-one POVM's are preprocessing clean. 

Notice that cleanness under preprocessing and extremality are properties com- 
pletely unrelated. Consider, e. g., the following rank-one POVM 

i|l)(l|,i|l)(l|,|2)(2|,...,|d)(d|. (4.4) 

The redundantly doubled outcome |1)(1| suggests at first sight that such a POVM 
cannot be extremal, namely, it cannot be the solution of any optimization problem. 
In this sense, such POVM "is not good". However, being rank-one, it is clean under 
both preprocessing and postprocessing. 

4.1.3 Positive maps 

Since now, we introduced two pseudo-orderings on the set of POVM's, the pre- 
processing ordering >~, and the postprocessing ordering >~p. In this Subsection, we 
will introduce two additional relations which can be established among POVM's, 
namely 

Definition 4.1.10 (positive preprocessing) We write P ^ Q and say that P 
is cleaner than Q under positive preprocessing, if and only if there exists a positive 
(non necessarily completely positive) map V such that Q = 'P(P). 

Definition 4.1.11 (range-inclusion) We write P Dj. Q and say that P range- 
includes Q, if and only z/Rng(Q) C Rng(P), where the range of a POVM is defined 
in Definition \1.3.1\ 

We simply have the following hierarchy of relations 

P^Q =^ P>Q =^ Pd, Q. (4.5) 

The converse is not always true. However, we have some results providing sufficient 
conditions for which some of the relations in Eq. ()4.5j) can be inverted. Proofs are 
very technical and can be found in Ref. [IH]. Here we just give the statements. 
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Theorem 4.1.12 Consider two POVM's P and Q with the same number of out- 
comes. Then the following statements are equivalent: 

1. P Dr Q 

2. There is a (unique) positive map £ : Span(P) — > Span(Q) with S(P) = Q. 

Notice tliat point (0) does not say tliat P ^ Q, since the positive map is defined 
only from Span(P) to Span(Q), and generally cannot be extended to a positive 
map on all B{J^). The following Theorems describe some situations in which it 
possible to extend the map £ to a positive map over all B(J^). 

Theorem 4.1.13 Consider two POVM's P and Q with the same number of out- 
comes. Then the following statements are equivalent: 

1. PyQ 

2. There is an informationally complete POVMlSA such that P(8)M Dr Q®M. 

3. P (g) M D,. Q O M holds for all POVM's M. 

Theorem 4.1.14 (abelian POVM) Consider two POVM's P and Q with the 
same number of outcomes. Let Q be abelian, namely QiQj = QjQi for all i,j. 
Then P Dr Q =^ P Q, and Eq. \4.5^ becomes a chain of equivalences. 

4.2 Inverting decoherence 

We will now focus our attention on a particularly nasty preprocessing of input 
states, namely on decoherence. Decoherence is universally considered, on one side, 
as the major practical limitation for communication and processing of quantum 
information. On the other side, decoherence yields the key concept to explain 
the transition from quantum to classical world [IH] due to the uncontrolled and 
unavoidable interactions with the environment. Great effort in the literature has 
been devoted to combat the effect of decoherence by engineering robust encoding- 
decoding schemes. Some authors have recently addressed a different approach to 
undo quantum noises by extracting classical information from the environment jSOj 
and exploiting it as an additional amount of side information useful to improve 
quantum communication performances |51j . 
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The recovery of quantum coherence from the environment is often a difficult 
task, e. g. when the environment is "too big" to be controlled, as for spontaneous 
emission of radiation. By regaining control on the environment the recovery can 
sometimes be actually accomplished, for example by keeping the emitted radiation 
inside a cavity. However, in some cases, the full recovery of quantum coherence 
becomes impossible even in principle, namely even when one has complete access 
to the environment. This naturally leads us to pose the following question: in 
which physical situations is possible to perfectly recover quantum coherence by 
monitoring the environment? 

4.2.1 Convex structure of decoherence maps 

A completely decohering evolution asymptotically cancels any quantum superpo- 
sition when reaching the stationary state, making any state diagonal in some fixed 
orthonormal basis — the basis depending on the particular system-environment in- 
teraction. In the Heisenberg picture we say that such a completely decohering 
evolution asymptotically maps the whole algebra of quantum observables into a 
"maximal classical algebra", that is a maximal set of commuting — namely jointly 
measureable — observables. Let's denote by 21^ the "quantum algebra" of all bounded 
operators B{Jif) on the finite dimensional Hilbert space J^, and by 2lc the "clas- 
sical algebra", namely any maximal Abelian subalgebra 2tc C 2lg. Clearly, all op- 
erators in 2lc can be jointly diagonalized on a common orthonormal basis, which 
in the following will be denoted as b = = l,...,d}. Then, the classical 

algebra 2lc is also the linear span of the one-dimensional projectors \k){k\, whence 
2lc is a d-dimensional vector space. According to the above general framework, 
we call (complete) decoherence map a completely positive identity-preserving (i. e. 
trace-preserving in the Schrodinger picture, see Subsection ll.4.2|) map S'^ which 
asymptotically maps any observable O G 2lg to a corresponding "classical observ- 
able" in 2lc, namely such that the limit lim„^oo(i^^)"(0) exists and belongs to the 
classical algebra 2tc for any O G 2lg. Here we denote with {S'^)^ the n-th iteration 
of the map S, implicitly assuming markovian evolution. 

It is easy to see that the set of decoherence maps is convex. The following 
Theorem shows that such maps enjoy a remarkably simple form: 
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Theorem 4.2.1 (Schur form) A map £^ preserves all elements of the maximal 
classical algebra 2lc if and only if it has the form 

£''{0)=ioO, (4.6) 

AoB denoting the Schur product of operators A andB, i. e. AoB = Xlfc /=i ^H-Bfe^lA;) (/|, 
{Afei} and {Bki} being the matrix elements of A and B in the basis h, and C,ki be- 
ing a correlation matrix, i.e. a positive semidefinite matrix with ^kk = 1 for all 
k = 1, . . . ,d. 

Theorem 14.2.11 states a hnear correspondence between maps preserving 2tc and 
correlation matrices, whence the two sets share the same convex structure. Then 
the map is extremal if and only if its correlation matrix is extremal. 

Since now we dealt with the dual map S'^ on bounded operators. The action 
of a decoherence map on quantum states is given in Schrodinger picture by 

£ip)=eop, (4.7) 

where T denotes transposition with respect to the basis b (also is a correlation 
matrix, hence in the following, we will drop the symbol T at the exponent). As 
a consequence, one has exponential decay of the off-diagonal elements of p, since 
I [£^"'(p)]fc/| = l^fc/T ■ \Pki\- In other words, any initial state p decays exponentially 
towards the completely decohered state 

Poo = ^pkk\k){k\. (4.8) 

k 

In Ref. 07j, it is proved the following 

Lemma 4.2.2 A map £ is an extremal decoherence map if and only if it is ex- 
tremal in the set of all maps. 

As a consequence of Lemma 14.2.21 the convex structure of decoherence maps can 
be obtained by application of the well known Choi Theorem jTH], which states that 
the canonical Kraus operators^ {^i}-, 1 < i < r, of every extremal map are such 
that their products {ElEj}^ 1 < < are linearly independent. A relevant 
consequence of this characterization is the following 

^For the definition of canonical Kraus decomposition, see Subsection 1 1.4 .21 
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Theorem 4.2.3 If S is an extremal decoherence map, then r < \fd. For qubits 
and qutrits any decoherence map is then random-unitary. 

This means that for qubits and qutrits extremal decoherence maps are unitary 
maps, since they admits a Kraus representation containing only one operator. 
Hence, for qubits and qutrits, every decoherence map can be written as 

£ip)=J2P^U^pUj, (4.9) 

i 

for some commuting unitary operators Ui G 2tc and probability distribution pj. 

4.2.2 Correcting decoherence by measuring the environ- 
ment 

In Ref jSU] it is shown that the only channels that can be perfectly inverted 
by monitoring the environment are the random-unitary ones. Therefore, it follows 
that one can perfectly correct any decoherence map for qubits and qutrits by mon- 
itoring the environment. The correction is achieved by retrieving the index i in 
Eq. ()4.9p via a measurement on the environment, and then by applying the inverse 
of the unitary transformation Ui on the system. Therefore, the random-unitary 
map simply leaks H{pi) bits of classical information into the environment {H de- 
noting the Shannon entropy), and the effects of decoherence can be completely 
eliminated by recovering such classical information, without any prior knowledge 
about the input state. The fact that decoherence maps are necessarily random- 
unitary is true only for qubits and qutrits. A counterexample in dimension d = 4 
can be found in Ref. jlZj. Such extremal decoherence maps with r > 2 represent 
a process which is fundamentally different from the random unitary one, corre- 
sponding to a leak of quantum information from the system to the environment, 
information that cannot be perfectly recovered from the environment [50] . 

Now we address the problem of estimating the amount of classical information 
needed in order to invert a random-unitary decoherence map. If the environment 
is initially in a pure state, say |0)e, a useful quantity to deal with is the so-called 
entropy exchange |S2] 5'ex defined as 



(4.10) 
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where is tlie reduced environment state after the interaction with the system 
in the state p, and S{p) = — Tr[plogp] is the von Neumann entropy. In the case 
of initially pure environment, the entropy exchange depends only on the map 
£ and on the input state of the system p, regardless of the particular system- 
environment interaction chosen to model S. It quantifies the information flow from 
the system to the environment and, for all input states p, one has the bound [221 
\S{S{p)) — S{p)\ < Sex{p), namely the entropy exchange Sex bounds the entropy 
production at each step of the decoherence process. 

In order to explicitly evaluate the entropy exchange for a decoherence process, 
we can then exploit a particular model interaction between system and environ- 
ment. This can be done noticing that it is always possible to write = {Gi\ek) 
for a suitable set of normalized vectors {|efc)}. Then, the map S{p) = ^ o p can 
be realized as S{p) = Tre[f/(p (g) |0) (0|e)?7^], where the unitary interaction U gives 
the transformation 



whence the final reduced state of the environment is = J2kPkk\^k){Gk\- Then, 
in order to evaluate S^x for a decoherence map S{p) = ^ op, it is possible to bypass 
the evaluation of the states \ei) of the environment, using the formula 



which follows immediately from the fact that y/Poo^y/poo, and are both reduced 
states of the same bipartite pure state y/pii\i)\Gi)- 




Notice that the unitary interaction U in Eq. ()4.1ip generalizes the usual form 
considered for quantum measurements with the quantum system interacting 
with a pointer, which is left in one of the (nonorthogonal) states {|efc)}. The more 
the pointer states are "classical" — i. e. distinguishable — the larger is the entropy 
exchange, whence the faster is the decoherence process. In the limit of orthogonal 
states, decoherence is istantaneous, i. e. S{p) = p^o- 

When a map can be inverted by monitoring the environment — i. e. in the 
random-unitary case — the entropy exchange Sex{I/d) provides a lower bound to 
the amount of classical information that must be collected from the environ- 
ment in order to perform the correction scheme of Ref. jSO]- In fact, assum- 
ing a random- unitary decomposition ()4.9p and using the formula Sex{p) = 



U\k) ® |0)e = \k) ® |efc). 



(4.11) 




(4.12) 
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S (Eij VPiPjT^AUipU^j]\i){j\), we obtain 

ScM/d) < H{pi). (4.13) 

The inequahty comes from the fact that the diagonal entries of a density matrix 
are always majorized by its eigenvalues jHH|, and it becomes equality if and only 
if Tr[f/i[/j]/d = i. e. the map admits a random- unitary decomposition with 
orthogonal unitary operators. Moreover, from Eq. ()4.12|1 we have Sex{I/d) = 

In Ref. |17j, it is proved that, for qubits, 5'(,^/2) quantifies exactly the minimum 
amount of classical information which must be extracted from the environment, 
while, for dimension d > 2, the bound in Eq. ()4.13|1 is generally strict and a 
counterexample is given for dimension d = 3. Notice that the same decoherence 
map may be obtainable by different random-unitary decompositions with different 
probability distributions {pi}, corresponding to different values of the information 
H[pi). However, for qubits it is always possible to perform a suitable measure- 
ment on the environment and to invert the decoherence map retrieving the minimal 
amount of information from the environment, namely 5'(.^/2). For example, con- 
sider the so-called random phase-kick model jHHl for decoherence of qubits 

1 r+oo 

Sip) = / e^^'^^/Ve-''"^/'e-^'/^M^, (4.14) 

V47rA J-oo 

which can be rewritten as 

£ip) = ^op, e= (^J:, (4.15) 

From Eq. ()4.14|1 . one could infer that the amount S of classical information that 
must be extracted from the environment is equal to the differential entropy of the 
gaussian probability density according to which the system is random phase-kicked, 
namely (see Ref. |31] ) 

growing logarithmically with A. This is actually not correct, since the minimum 
amount of classical information needed is 

1 - 

S" = -plogsp- (1 -p)log2(l -p), p = . (4.17) 
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Figure 4.2: Tlie amount of classical information 5*, expressed in bits, leaking into 
the environment at every application of the random phase-kick model of decoher- 
ence for qubits, as function of the parameter A, see Eqs. fl4.14|l and ()4.15j) . S 
tends to the limit value of 1 bit, since every qubit decoherence map can be written 
as a random-unitary process involving only two unitaries (see the footnote in the 
previous page). 



l+e- 



In fact, ^/2 in Eq. ()4.15p can be simply diagonalized and has eigenvalues | , ^ 
In Figure I3I21 there is the plot of the amount of classical information S in Eq. ()4.17|) 
as a function of the parameter A modelling the decoherence rate in Eqs. ()4.14|1 
and ()4.15j) . The curve tends to the finite limit of one bit, contrarily to what 
happens in Eq. ()4.1fij) . 
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